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CHAPTER I 

1.1 ADMINISTRATIVE INFORMATION 

This report covers the work accomplished on N.A.S.A, 
contract number NAS8-31233 during the time period 11/1/74 
through 1/31/76 for tnc. Oeorge C. Marshall Space Flight 
Center, Alabama 35812. The principal investigator for 
this study was Dr. Dara W, Childs, an Associate Professor 
of Mechanical Engineering at the University of Louisville. 
Direction for the study has been provided by Mr. Loren 
Gross and Mr. Woody Pitsenberger , Marshall Space Flight 
Center. 

1.2 RESEARCH OBJECTIVES AND ABSTPvACT 

The rotordy naitii c behavior of the SSME (Space Shuttle 
Main Engine) KPOTP (High Pressure Oxygen Turbopump) and 
HPFTP (High Pressure Fuel Pump) have been analyzed in this 
study, with the follovjing objectives: 

(a) The identification of potential rotordynamic problem 
areas vjhich might arise during operation of these 
units prior to their testing. 

(b) The investigation of alternative procedures for cor- 
recting potential rotordynamic problems should they 
occur . 

(c) The development of an adequate analytic and physical 
understanding of the turbopump rotordynamics to im- 
prove the probability of a correct diagnosis of ro- 
tordynamic problems from test data. 
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Improved transient rotordynandc models are developed 
in this study for both turbopumps based on 111** The pre- 
sent study is a continuation of prior rotordynamic studies 
for the SSME turbopumps [2^3], and is concerned primarily 
with rotordynamic stability. The transient models devel- 
oped here employ the results of Black [4-7] and Yamada [81 
to model the hydrodynamic forces developed by the SSME tur- 
bopump seals. Linear stability analysis is performed for 
the SSME turbop’omp rotordynamlcs models, which includes 
gyroscopic effects , seal forces , speed-dependent bearing 
characteristics, and internal rotor damping. 

The stability results for the HPPTP are disturbing in 
that the first mode of the rotor is predicted to be unsta- 
ble due to seal forces alone with zero internal rotor dam- 
ping. The addition of internal rotor damping on the order 
of one percent for the first two modes significantly wor- 
sens the situation. 


Identifies entry in Reference Section 
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The following two i*sdesigns were considered for the 
puinp and turbine-end bearing carrier structures to control 
the predicted rotordynamic instability: 

(a) Stiffness asymmetry in the bearing carrier structure 
was examined/ and significantly improves the stability 
characteristics of the rotor. Specifically/ a bearing 
carrier would be redesigned so that its radial stiff- 
ness in two axial orthogonal planes are markedly dif- 
ferent. The carriers would then be assembled in the 
turbopumps so that their major and minor stiffness 
planes coincide/ preferably with the major stiffness 
plane aligned with the hydrodynamic sideload to reduce 
steady-state rotor deflections. The adequacy of this 
approach depends on the degree of rotor instability as 
determined by hardware testing. 

(b) A redesign was also considered for the bearing carrier 
structures to develops squeeze-film damper surfaces. 

The results of this phase of the study are mixed. The 
squeeze-film-dampers improve rotor stability/ and in 
some cases are sufficient to contain or eliminate the 
rotor instability. However, their capacity in this re- 
gard is limited/ and they can be overloaded. Kirk and 
Gunter's nonlinear model for squeeze-film dampers 1173 
was used to model the squeeze-film damper forces. 

A combination of the two redesign approaches was also 
considered, again with mixed results. The adequacy of a 


4 
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combined approach depends on the degree of rotor instabil- 
ity, and the net improvement from a combination of the two 
approaches appears to be less than the sum of their im- 
provements taken separately. It is clearly possible to 
overload a combined stiffness asymmetry/squeeze-film dam- 
per design. 

Stability results for the HPOTP rotor indicates that 
it is more stable than the HPFTP rotor, but is also poten- 
tially subject to sub synchronous whirling. Transient si- 
mulation models have been developed for the analysis of 
unstable motion of this rotor if test results of the rotor 
confirm the presence of an unstable subsynchronous whirling 
mode. 
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CHAPTER II 

ROTORDYNAMIC SEAL ANALYSIS FOR SSME TURBOPUMPS 
2.1 THE ANALYSIS OF YAMADA AND BLACK ET AL. 


In this study, the work of Black et al [4-7] was used 
to define the load-deflection properties of axial-flow 
sealing surfaces in both turbopumps. Black’s analysis is 
in turn based on a seal-flov7 analysis by Yamada [8] , devel- 
oped for a "plain" seal, resembling a journal bearing in 
the sense that the seal has neither steps nor annular ser- 
rations. A reviev; of the Yamada-Black procedure follows. 

The implementation of Black's analysis for load-de- 
flection properties of seals requires the definition of 
their flov/ characteristics. Black employs the following 
relationship from Yamada [8] . 


AP = (1 + C + '2a)pv^/2 = Cctpv2/2 (1) 

where C is a (constant) entry- loss coefficient, p is the 
fluid density, v is the average fluid velocity, and a is 
a friction-loss coefficient defined by 

a = \l/& (2) 

In the above, Jl is the seal length, 6 is the radial clear- 
ance, and X has been defined by Yamada to be the following 
function of the axial and radial Reynold's numbers (R^ , R„) 

a, JT 


X = 0,079 R. 


;?A[i - 


(3> 


R^ = 2v6/v, R^ = r(^6/v 


v/here v is the fluid's kinematic viscosity, and r is the 
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seal radius. 


Black's analysis yields a load-deflection definition 
of the form 



where | is the rotor spin velocity, and 


K = (5) 

C = lil®? c = P 2 I r m 
T = A/V 

In addition, the coefficients ]iq , m, y 2 , are a function 
of the coefficients C, a, and $, where 

a = (J|21)2/U + (6) 

The factor 6 is included in the analysis to account for the 
change in the friction factor X due to a perturbed displace- 
ment of the rotor away from a centered position in the seal. 
Plots of yo / 111/ V2 s^^re provided in Fig. 1 for 5 = 0.5 
(Ref. [5] ) . These coefficients are comparatively insensi- 
tive to reasonable variations in 


i 
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Black's analysis accounts for the length of the seal 
via the following formulae 
Po (e) = po/(l + 0.28 e^) 

yi(e) = Pi/(1 + 0.23 e2) (7) 

V2 0.06 e^) 

where e = Jl/r. 

2.2 SSME SEAL-FLOW ANALYSIS PROCEDURE 

As noted previously/ Yamada's flow analysis was devel- 
oped for "plain" seals / which basically look like journal 
bearings. The seals in the high pressure SSME turbopximps 
are not generally of this configuration, since they employ 
both steps and annular serrations. This section deals with 
the analysis procedure which was developed to "fit" the re- 
sults of Yamada's analysis to the SSME seals. 

Yamada's analysis can be used to iteratively calculate 
the leakage flowrate through a seal for a specified differ- 
ential pressure AP and a specified entry loss coefficient 
C. Specifically, one assumes a value for X, calculates a 
from Eq. (2) , v from Eq, (1) , and then calculates a nev7 
value for X from Eq. (3) . This procedure can be repeated 
until "stable" values for X and v are achieved- However, 
this procedure was not initially followed in the present 
study. Instead, Rocketdyne personnel supplied leakage 
flowrates, and differential pressure data for each seal at 
FPL. The calculations which yielded this data were based 
on Rocketdyne' s considerable past experiences with the 
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type of seals used in the turbopumps , and accounted for the 
reduction in radial clearances due to the centrifugally in- 
duced growth of the rotating assembly. 

The assumption was made in fitting the Rocketdyne FPL 
leakage data to the Yamada seal model that the Yaraada re- 
sults were valid for the constant-clearance, constant-radi- 
us segments of seals- The data -were fitted to the ‘ - ; • 

model by calculating the (average) velocity at each segment 
v;ithin the seal, and distributing the total iP across the 
seal to each segment on a v^ basis. ■■ Then for each segment within 
the seal, Eq. (3) is used to calculate a, and the entrance 
loss coefficient ? is calculated from Eq. (1) ' to satisfy 
AP and v^ (Rocketdyne data) and a (Yamada seal model) . 

The clearances of the seals at FPL speeds used in this 
study coincided with estimates supplied by Rocketdyne. 

Where estimates were not provided, the clearances v/ere simply 
adjusted to yield "reasonable” entry loss coefficients. This 
adjustment is fairly simple, since for a specified flowrate, 
the axial Reynolds number R is independent of the clear- 
ance 6. Thus, X (and a) is also largely independent of 6 , 
and Eq. (1) can be used to calculate a suitable average 
flow velocity v, and 6 is selected to yield the prescribed 
flow rate. 

The above procedures yield entrance loss coefficients 
and friction factors for each segment of each seal at FPL. 

The follov7ing procedure was followed in extrapolating this 
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data across the major portion of the operating speed range 
of the turbopumps. The following relationships were as- 
sumed for the differential pressure across a seal and the 
seals radial clearance 

&P = d^2, fr = «o “ (8) 

where Sq the zero-speed clearance. The constants d and 
b are calculated for each seal from the known FPL condi- 
tions. hence ^ the values of AP and S are readily calculated 
for any specified running speed. The entrance-loss coef- 
ficients which were calculated for FPL speeds were used at 
all speeds, and the flow-rate for a specified was calcu- 
lated in the iterative fashion outlined above. Specifi- 
cally, for a given segment within the seal, AP and 6 V7ere 
calculated from Eq. (8) , an initial value of A was assigned, 
and an average velocity calculated from Eq. (1) . Eq. (3) 
is then used to calculate a new X, and the cycle is repeated 
until consistent values of X and v are obtained. The flow- 
rate is then calculated from v and a . 

The viscosity and density of the fluid within the seal 
are required for Yamada’s seal analysis. The required data 
for hydrogen and oxygen were taken from References [9] and 
[101 / respectj-vely , for pressure and temperature data sup- 
plied by Rocketdyne personnel. Fluid properties in the tur- 
bine areas of the turbopumps were taken from calcu] ations by 
K. Gross of Marshall Space Flight Center. 

The assumptions used here in defining the flow proper- 
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ties of seals are generally supported by the work of Isaac- 
son as reported by Fasheh [11] . Isaacson tested a range of 
seal teeth configuration and reached the following conclu- 
sions: 

(a) stepped type labyrinth seals are preferred for 
leakage control, 

(b) stepped type seals are comparatively insensj.tive 
to changes in clearances and eccentricities, and 

(c) rotational effects on leakage are reduced for 1 
large AP. 

A stepped lab;* rinth seal is designed to achieve f lov7 con- 
striction primarily by changes in the direction of the flow 
as opposed to fluid friction losses arising due to viscosity. 
Hence, conclusion (b) would be anticipated for a properly 
designed seal, and the implication for the seal model of Eq. 
(1) is that the friction loss term 2a would be small in com- 
parison to(l+g). Further, Isaacson's test results indicate 
that the over all seal coefficient (for stepped labyrinth 
seals) decreases as the pjxial Reynolds number is in- 
creased, and tends to approach an asymptotic value for 
on the order of .5 x 10^. This characteristic is interpre- 
ted here to mean that the entrance loss coefficient ^ is in- 
dependent of R^, and the reduction in results from chan- 
ges in X (and hence a) as defined by Eq . (3). 


12 


2.3 SEAL ANALYSIS FOR THE SSME HPFTP 

The analysis outlined in the preceeding section was car- 
ried out for the HPFTP seals illustrated in figure 2, for 
the following speed set; 

= 1654 rad/sec {15 ,.800 rpm) 

I 2 = 2107 rad/sec (20,100' rpm) 

^3 = 2560 rad/sec (24,500 rpm) 

( 9 ) 

= 3016 rad/sec (28,800 rpm) 

Is - 3469 rad/sec (33,200 rpm) 

Is = 3920 rad/sec (37,470 rpm) = FPL 

The speed range v;as selected to include the speed at FPL 

(le = 3920 rad/sec) and the second critical speed. Table 1 
contains a representative set of data and results at FPL. 

The entries in the second column of this table is the num- 
ber of constant-clearance constant-radius annular segments used 
in modeling each seal. The entries in the AP column were 
supplied by Rocketdyne, and the radial clearances at FPL are 
generally based on Rocketdyne estimates. The axial position 
entries in Table 1 are with respect to the pump end of the 
turbopump . 

The viscosity, specific v^eight, and flowrate of the 
fluid VJithin the seals for the speed set of Eq. (9) are gi- 
ven in Tables 2, 3, and 4, respectively. The FPL flowrate 
in Table 4 \jas provided by Rocketdyne. The remaining flow- 
rates vjere calculated by the method outlined in the prece- 
ding section. The results for some of these seals are con- 



TABLE 1: HPFTP Seal Data and Results at FPL 


Seal No- Axial 


No. 

Segments 

Position (in. ) 

AP (Psi) 

6 (in. ) 

S 

a 

R_ X 10 ® 

R_ X 10' 

1 

4 

5.25 

1325. 

.005 

0.51 

0.025 

2.01 

2.21 

2 

3 

7.42 

980. 

.005 

0.99 

0.029 

3.29 

2.58 

3 

3 

8.50 

1980. 

.005 

0.38 

0.265 

1.27 

4.28 

4 

4 

10 . 16 

1153. 

.005 

0.32 

0.025 

2.67 

2.89 

D 

3 

12.37 . 

956. 

.005 

0.95 

0.030 

3.10 

2.34 

6 

3 

13.33 

1990. 

.005 

0.53 

0.315 

1.14 

3.54 

7 

4 

15-13 

1319. 

.005 

0.43 

0.024 

2.38 

2.64 

8 

9 

20.37 

335. 

.005 

1.48 

0.014 

1.48 

0.908 

9 

1 

22.25 

359. 

.010 

0.39 

0.045 

0.671 

1.26 

10 

1 

24.88 

485. 

.005 

0,67 

0.062 

0.702 

1.25 

11 

3 

25.79 

245- 

.006 

0.28 

0.125 

1.09 

2.42 


I 
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Table 2; HPFTP Seal Viscosity (Ib'n/ft sec) x 10® at the 
Selected Operating Speeds of Eq. lyj 


Seal 

No. 

FPL 

+5 

I4 

I3 



1 

1.150 

1.088 

1.029 

0.969 

0.924 

0,887 

2 

0.827 

0.760 

0.713 

.0.657 

0.611 

0.570 

3 

0.710 

0.644 

0.578 

0.517 

0.448 

0.398 

4 

0.890 

0.790 

0.682 

0.548 

0.483 

0.356 

5 

0.330 

0.816 

0-687 

0.548 

0.499 

0.384 

6 

0.852 

0.740 

0.636 

0.530 

0.426 

0.310 

7 

1.002 

0.851 

0.720 

0.594 

0.478 

0.352 

B 

0.698 

0.616 

0.576 

0.410- 

0.341 

0.275 

9 

1.750 

1.570 

1.450 

1.34 

1.23 

1.10 

10 

1.780 

1.630 

1.500 

1.38 

1.26 

1.11 

11 

0.550 

0.487 

0.427 

0.376 

0.335 

0.305 


tPable 3: HPFTP, Seal Specific VJeight (Ib/ft^} at the 

Operating Speeds of Eq* (9) 




Seal 

No. 

FPL 

h 

• 

<{ii, 

^3 

^2 



1 

4.74 

4.67 

4.61 

4.54 

4.48 

4.44 


2 

4.45 

4.34 

4.23 

4.10 

3.98 

3.86 


3 

4.25 

4.09 

3.90 

3.69 

3.41 

2.87 


4 

4.90 

4.50 

4.20 

3.79 

3.54 

2.92 

• 

5 

4.73 

4.51 

4.21 

3.97 

3.61 

3.04 

* 

6 

4.60 

4.30 

4.05 

3.69 

3.21 

2.40 


7 

4.90 

4.60 

4.31 

3.90 

3.42 

2.65 

h 

8 

4,20 

3.86 

3.74 

2.88 

2.22 

1.94 

\ 

9 

0.78 

0.686 

0.550 

0.411 

0.294 

0.138 


10 

0.99 

0.833 

0.664 

0.493 

0.348 

0.154 


11 

3.40 

3.01 

2.530 

2.02 

1.45 

0.927 

- 
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Table 4: HPPTP Flow Rate (Ib/sec) at the Operating 

Speeds of Eg. (9) 


Seal 


• 

■ 

* 

• 


No. 

FPL 

'i>S 

4't* 

‘l'3, 

^2- 

ii. 

1 

2.383 

2.8794 

3.0853 

3.0346 

2.7769 

2.3561 

2 

2.519 

4.6292 

5.8077 

6.1796 

5.9050 

5.1217 

3 

1.239 

1.2665 

1.2187 

1.1066 

0.9415 

0.7164 

4 

2.404 

3.3473 

3.7796 

3.7230 

3.4076 

2.6860 

5 

2.587 

4.7021 

5.7743 

6.0609 

5.6055 

4.5309 

6 

1.229 

1.2380 

1.1833 

1.0559 

0.8721 

0.6254 

7 

2.473 

3.4805 

3.9368 

3.8827 

3.4435 

2.6306 

8 

0.301 

0-26096 

0.22753 

0.17224 

0.1261 

0.0935 

9 

3.200 

2.8640 

2.3712 

1. 8289 

1.3238 

0.7326 

10 

3.200 

3.3648 

3.1308 

2.6173 

1.9972 

1.1204 

11 

0.375 

0.31198 

0.24830 

0.18779 

0.1302 

0.0810 
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trary to one's normal expectations, since the flowrate in- 
creases as speed and AP decrease. However, for these seals, 
the radial clearance is simply increasing faster than AP is 
uecreasing , with a resultant increase in flow. 

The radii and lengths of seal segments were taken from 
component drawings. The lengths of the segments within the 
seal are summarized below: 

SEAL LENGTHS (in.) 


1 

^i 

= 

.'030; 

i = 1, 2, 3, 4 


2 

^i 

= 

.030; 

i == 1, 2, 3 


3 


= 

.420, 

Z 2 ~ • 470 , 3 

= .516 

4 

£. 

r 

- 

.030; 

i=l, 2, 3, 4 


5 

1 . 

X 

- 

.030; 

i = 1, 2, 3 


6 

S-i 


.420, 

Z 2 ~ • 470 , Z 2 

= .516 

7 

z. 

X 

= 

.030; 

i = 1 , 2 , 3, 4 


8 

z. 

X 

— 

. 100 ; 

i = 1 , 2 , . . 

.9 

9 


= 

.099 



10 

«■! 


.1154 



11 

^1 

-- 

.215, 

Z 2 = .200 , A 3 

= .135 


The segments in seals 3 and 6 are seen to be significantly 

longer than the segments of the other seals. The segments | 

of these seals closely resemble seals tested by Black and | 

Cochrane [7] in that each segment is finely grooved with com- ] 

! paratively shallow annular serrations. Black and Cochrane , 

suggest for this type of seal that the total length of the 
seal (including serrations) be used. The other seals of the 
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turbopump generally differ from seals 3 and 6 in that the 
axial flow segments are separated by either changes in rad- 
ius or comparatively deep serrations. Hence, for the re- 
maining seals, only the lengths of the segments themselves 
(excluding serrations) were included in the analysis. 

Returning to Table 1, the last four columns contain 
average values for a, R^, and at FPL for the segments 
within seals. The values for S can be compared to a range 
of 0.1 to 0.5 cited by Black. The Reynolds numbers in this 
chart are considerably greater than those achieved by Black 
experimentally. The a values for seals 3 and 6 are approx- 
imately an order of magnitude higher than any other seals, 
because of the comparatively long lengths of the segments 
v^ithin these seals. 

Tables 5 through 8 contain the physical stiffness and 

damping coefficients for the HPPTP seals for the speed set of 

Eg. (9) . Prom Eq. (4) these coefficients are defined by 

K = RCirrAP/A), k = MirrAP/l) 

- ^ - CIO) 

C = C(iTriP/A), c = c(urAP/A) 

In calculating these coefficients, the seal segments were 
assumed to be in parallel, i.e., stiffness and damping coef- 
ficients were calculated for individual segments of the 
seal, and then summed to obtain effective coefficients for 
the seal. An inspection of the results of Tables 5 through 
8 demon? -ates that the diagonal coefficients are approxi- 
mately ct'. order of magnitude greater than the off diagonal 


1 

1 
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Table 5: HPFTP Seal Diagonal Stiffness Coefficient 

K(lb/in) at the Operating Speeds of Eq, (9) 


Seal 

• 

• 

• 

• 

« 

No. 

FPL 


4’4 


4>2 


1 

1832. 

731. 

352. 

185. 

101. 

54.5 

2 

1417. 

219. 

72.4 

31.0 

14.9 

7.42 

3 

1.223x105 

7-658x10*' 4 

.840x10" 

3.035x10*^ 

1.851x10" 

1.053x10" 

4 

169B. 

454. 

179. 

81.5 

41.1 

20.0 

5 

1460. 

223. 

71.4 

29.4 

14.1 

6-. 61 . 

6 

1.193x10^ 

7.402x10** 4 

.634x10" 

2.872x10'' 

1.734x10" 

9825. 

7 

1826- 

483. 

189. 

86.. 8 

42.8 

20.9 

Q 

26.16 . 

19.3 

.14.2 

- 9.20 

5.95 

in 

fn 

3 

2218. 

1501. 

1025. 

692. 

454. 

290- 

10 

4621. 

2173. 

1153. 

652. 

377. 

221. 

11 

2183. 

1719, 

1313, 

963. 

674. 

438. 


Table 

6: HPFTP Off 

-Diagonal 

'Seal Stiffness Coefficient 



k(lb/in) 

at the Operating Speeds of Eq. 

(9) 

Seal 






No. 

FPL 

$5 

^4 

^3 

|'2 

la 

1 

275.4 

159. 

97.6 

60.1 

36.3 

20.5 

2 

168.2 

64.1' 

33.0 

18.5 

10.5 

5.67 

3 

1-870x10*' 

1.343x10" 

93.36. 

6206. 

3866. 

' 2111. 

4 

315.4 

150. 

82.6 

46.6 

26.6 

13.5 

5 

175,6 

66.3 

33.4 

18.4 * 

10.1 

5.11 

6 

1.747x10** 

1.230xl0" 

8480. 

5534. 

3346. 

' 1720. 

7 

300.9 

144. 

79.6 

44.9 

24,9 

12.2 

8 

5.492 

4.04 

2. .95 

1.04 

1.08 

0.615 

9 

261.8 

180. 

115. 

68,5 

37.6 

15.3 

10 ■ 

409.1 

237. 

136. 

74,9 

38,8 

14.8 

11 

859.7 

631. 

434. 

'■276. 

115, 

73.9 
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Table 7: HPFTP Diagonal Seal Damping Coefficient 

C(lb/in/sec) at the Operating Speeds of Eq. (9) 


Seal 


No. 

FPL 

h 

^4 

• 

<1>3 

• 

4>2 


1 

0.1^05 

0.919 

.0647 

.0470 

.0344 

.0248 

2 

0.0858 

.0370 

.0219 

0.0144 

0.00995 

0.00686 

3 

9.543 

7.74 

6.14 

4-85 

3.67 

2.55 

4 

0.1609 

0.0864 

0.0547 

0.0364 

0.0253 

0.0163 

5 

0.0896 

0.0382 

0.0222 

0.C144 

0.00961 

0.00618 

6 

8.915 

7.09 

5.62 

4.32 

3.18 

2.08 

7 

0.1535 

0.0816 

0.0528 

0.0351 

0.0236 

0.0148 

8 

0.0028 

0.00233 

0-00196 

0.00144 

0.00102 

0 .744x10“ 

9 

0.1336 

0.104 

0,0765 

0.0535 

0.0357 

0.0185 

10 

0.2087 

0.137 

0.0904 

0.0585 

0.0369 

0.0179 

11 

0.4386 

0.364 

0.288 

0.216 

0.147 

0.0894 




Table 8: HPFTP Off-Diagonal Seal Damping Coefficient 

c(lb/in/sec) at Selected Operating Speeds 


Seal 

No. 

FPL 

• 

*1*5 

* 

m 

^*3 

• 

<1*2 


i 

0.820x10“^* 

o.ssgxio"'^ 

0.724X10"** 

0.142x10"** 

0.962x10“** 

0.665x10“** 

2 

0.2862x10'*3 

0.l09xl0“3 

'0.618x10*''* 

0. 396x10"** 

0.267x10“** 

0.181x10“** 

3 

0.4876 

0.326 

0.223 

0.151 

0.105 

0.0633 

4 

-0.1405x10“3 

-0.907xl0“‘» 

-0.597x10“** 

-0.390x10“** 

-0.266x10"** 

-0.158x10"** 

5 

0.2965x10" 3 

0 . 109x10"^ 

0.593x10"** 

0-369x10“** 

0.223x10"** 

0.137x10"** 

6 

0.6232 

0.417 

0.289 

0.195 

0.126 

0.0688 

7 

0.3364x10"^ 

-0.163x10“*^ 

-0.139x10“** 

-O. 101x10"** 

-0.679x10-5 

-0.389x10"5 

8 

0.1404x10“*^ 

0.122x10“^* 

0.922x10“^ 

0.591x1C"5 

0.370x10"5 

0.251x10“5 

9 

0.1014x10“'^ 

0.565x10“5 

0.999x10"6 

0.117x10“6 

0.390x10"® 

0.800x10"® 

10 

0.7076x10“3 

0.368x10“^ 

0.799x10“3 

0.106x10"3 

0.547x10"** 

0,178x10“** 

11 

0.9461x10“3 

o.esixio-^ 

0.757x10"3 

0.673x10“3 

0.564x1O“3 

0.4 17x10“ 3 


to 

H* 


Table 9: HPOTP Seal Data and Results at FPL 


Seal 

No. 

No. 

Segments 

Axial 

Position (in) 

AP (Psi) 

1 

5 

0.58 

2099 . 

2 

3 

2.64 

6581. 

3 

1 

19.71 

90.' 

4 

1 

20.48 

1865. 

5 

2 

21.44 

156. 


3 

22.47 

581. 

7 

2 

23.64 

153. 



C 

a 

R^xlO ^ 

R xlO" 
a 

006 

0.720 

0.137 

0.725 

0.605 

005 

0.78 

0.185 

0.128 

1.27 

0025 

0.88 

0.963 

0.0026 

0,0456 

0025 

1.31 

0.591 

0.0432 

0.859 

005 

1.13 

0.055 

0.406 

0-323 

0125 

0.654 

0.046 

1.01 

2.25 

005 

1.73 

0.063 

0.517 

0.317 


to 

to 
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coefficients/ i.e. , K > k; C > c. Further, the coeffi- 
cients for seals 3 and 6 are an approximate order of mag- 
nitude greater than the coefficients for the other seals. 

JPhis result is accounted for by the comparatively large val- 
ues of a for these seals . The negative values for the cross- 
coupling damping tejnns c in Table 8 arise because (contrary 
to the 1 ‘esults of figure 1 from Ref. [5] ) , for sufficiently 
small values of a, \i 2 is negative. However, the magnitude 
of these coefficients are so small that they have a. negligible 
effect oh. turbopump. rotordynam.ics irrespective of their sign. 

The seal inertia coefficients ra of Eq . (4) were not 

used in the analyses of this investigation. Preliminary 
calculations indicated that these coefficients were negli- 
gibly small for the HPFTP. 


2.4 SEAL ANALYSIS FOR THE SSME HPOTP 

The analysis outlined in Section 2.2 and applied in the pre- 
ceding section to the HPFTP v/as carried out for the HPOTP 
seals illustrated in figure 3 for the following speed set: 

= 1345 rad/sec = (12,860 rpra) 

<^2 = 1728 rad/sec = (16,400 rpm) 

(^3 = 2112 rad/sec = (20,200 rpm) 

( 11 ) 

= 2496 rad/sec = (23,900 rpm) 

^5 = 2879 rad/sec = (27,500 rpm) 

<f>6 = 3263 rad/sec = (31,200 rpm) = FPL 
The speed was selected to include the FPL speed (^g = 3263 
rad/sec) and the first critical speed. The axial positions 
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1 . 

of the seals with respect to the pump end, and the niamber : 

of segments used in modeling each seal are given in Table 9 . I 

The lengths of the segments within the seals are suiratiar- | 

ized below; ! 

f, 

SEAL LENGTHS (in-) | 

1 £. = .124; i = 1, ...5 1 ; 

I 

2 = .183, £2 = .229, £3 = .229 i 

■ ■ 

3 jll = .250 

4 £1 = .320 

5 2 - 1 “ .035, £2 ~ .035 

i 

6 ^i ~ .150; i = 1, 2, 3 ' 

7 “ .035, £2 = .035 I 

i ; 

The viscosity and specific weight of fluid v/ithin the se^ls ’ 

are given in Tables 10 and 11. The fluid within seals 1 and j | 

2 is lox, within seals 3 and 4 is hydrogen, and within seals | 

5 through 7 is a hot gas mixture of the combustion products . | | 

Seal flowrates at the speed set of Eq. (11) are provi- ; 

t M 

ded in Table 12, The flowrates at FPL in this chart are j 

Rocketdyne estimates. The remaining flov/rates were calcu- i 

lated by the analysis procedure of section 2.2 based on the I 

FPL differential pressure and radial clearance data of Table ^ 1 

9. ; U 

Table 9 also provides average values for C/ o , and j 

R^ for the seal segments at FPL. This data generally resem- ! 

! 

bles the data of Table 1 for the HPFTP, except for seal 3. j 

This seal has a low flov;rate, which yields low Reynold's ■ 

numbers, and a large X, The combination of a large X and a large 
length yields a large a . Seal 4 is the principal turbine 
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Table 10; HPOTP Seal Viscosity (Tbm/f t . sec . ) x 10 ^ 
at the Operating Speeds of Eq. (11) 


Seal 







No. 

FPL 

^5 


^ 3 

I 2 

^ 1 

1 

20. 

18.72 

16.60 

15.50 

14.50 

13.70 

2 

17. 

16.15 

15.18 

14.46 

13.81 

13.35 

3 

0.661 

0.684 

0.684 

0.684 

0.684 

0.684 

4 

0.683 

0.695 

0.695 

0.695 

0.690 

0.690 

5 

1.394 

1.304 

1.245 

1.187 

1.128 

1.039 

6 

1,09 

1.053 

1.037 

1.022 

1.006 

0.990 

7 

1.45 

1.341 

1.274 

1.207 

1.141 

1.074 


Table 11: HPOTP Seal Specific Weight (Ib/ft^) at the 



' 

Operating 

Speeds of 

Eq. (11) 



Seal 







No. 

FPL 

^5 


^3 



1 

75 : 

74.7 

74.5 

72.8 

7 2,2 

71.7 

2 

73. 

73.3 

72.7 

72.2 

. 7 

71.3 

3 

0.0238 

0.0113 

0.0078 

.0056 

0.0033 

0.0023 

4 

0.4006 

0.324 

0.240 

0.176 

0.122 

0.073 

5 

0.909 

0.710 

0.573 

0.435 

0.298 

0.160 

6 

1.358 

0.993 

0.770 

0.54 8 

0.325 

0.102 

7 

1.218 

0.861 

0.681 

0.500 

0.319 

0.139 


1 
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Table 12: KPOTP Seal Flowrate (Ib/sec) at the 

Operating Speeds -of Eg. (11) 


Seal 


No. 

PPL 

<^5 

(fij 

^3 

<{>2 

1 

4.5 

4.278 

3.942 

3.466 

2.936 

2 

9.5 

10.28 

10.21 

9.60 

8.492 

3 

0.01 

.006 

. 0043 

0.0031 

0.0019 

4 

0.20 

.160 

0.120 

0.087 

0.059 

5 

0.65 

0.657 

0.627 

0.534 

0.402 

6 

1.851 

1.360 

1.061 

0.771 

0.492 

7 

0.650 

0.680 

0.640 

0.537 

0.390 


125 

396 

0011 

031 

228 

195 

199 


Table 13: HPOTP Seal Diagonal Stiffness Coefficient 

K{lb/in) at 'the Operating Speeds of Eq; (11) 


Seal 


No. 

FPL 


• 

•J'tf 

* 

m 

$2 

<i» a 

1 

1.489x10** 

1,016x10** 

. 6810. 

4514. 

2875. 

1507. 

2 

1.418x105 

9.176x10** 

5.998x10** 

3.893x10** 

2.441x10** 

1.268x10** 

3 

5921. 

4524. 

3334. 

2343. 

1529. 

943.4 

4 

1.403x105 . 

1.089x105 

8.204x10** 

5.919x10** 

4.008x10** 

2.286x10** 

5 

564.0 

224.3 

110.5 

59.05 

32.48 

15.49 

6 

1973. 

1394. 

1049. 

769.7 

546.6 

324.1 

7 

487.1 

210.0 

102.7 

54.35 

29.67 

14.29 


'j 
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Table 14: HPOTP Seal Off-Diagonal Stiffness Coefficient 

k(lb/in) at the Operating Speeds of Eg. (11) 


Seal 


No. 

FPL 

4 

• 

m 

• 

^2 

^ 1 

1 

4968, 

3653. 

2620. 

1731. 

1147. 

509.2 

2 

2.462xl0't 

1.946x10“ 

1.454x10“ 

1.029x10“ 

6775. 

3187. 

3 

31.54 

15.69 

9.321 

5.392 

2.584 

1.058 

4 

1090. 

7391. 

461.4 

272.8 

145.9 

52.65 

5 

30.03 

14.41 

8.060 

4.390 

2.195 

0.705 

6 

474.6 

264.2 

170.5 

100.7 

50.72 

12.49 

7 

23.88 

13.79 

7.632 

4.085 

1.972 

6.570 



Table 15 : 

HPOTP Seal Diagonal 
C (Ib/i-n/sec) at the 

Damping Coefficient 
Operating Speeds of Eg . 

(11) 

Seal 

No. 

FPL 

•is 

• 

(j»4 

^3 

* 

4’2 

n 

1 

3.045 

2.538 

2.101 

1.687 

1.328 

0.757 

2 

15.09 

13.51 

11.66 

9.742 

7.841 

4.739 

3 

0.0193 

0.0109 

0.00747 

0.00510 

0.0030 

0.00157 

4 

0.668 

0.513 

0.370 

0.258 

0.168 

0.0783 

5 

0.0184 

0.0100 

0.00646 

0.00415 

0.00254 

0.00104 

6 

0.291 

0.183 

0.137 

0.0953 

0.0587 

0.0186 

7 

0.0147 

0.00957 

0.00611 

0.00387 

0.00228 

0.848x1 


Table 16; HPOTP Seal Off-Diagonal Damping 
Coefficient c{lb/in/sec) at the 
Operating Speeds of Eg. (11) 


Seal 

No. 

PPL 

^5 

K 

*3 

^2 



0.152 

0.120 

0.0953 

0.0740 

0.0574 

0.0574 

2 

0.803 

0.519 

0-333 

0.216 

0.140 

0,167 

3 

0.287x10"3 

0.118x10“3 

0.693x10“** 

0.416x10"** 

0.196x10“** 

0.111x10“'* 

4 

0.0105 

0.00739 

0.00470 

0.00290 

O. 00164 

0.769x10“** 

5 

0.753x10“'* 

0.390x10“** 

0.218x10“** 

0.120x10"** 

0.601x10“5 

0.251xl0“3 

6 

0.00164 

0-00142 

0.940x10“3 

0.563x10-3 

0.275x10-3 

0.864x10"** 

7 

0.111xl0“3 

0.486x10“** 

0.267x10"** 

0.142x10“** 

0.664x10“3 

0.216x10“3 


NJ 

VO 
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f lOating-ring seal, and at PPL is choked. Hence, the AP ci- 
ted in Table 9 accounts only for the pressure drop due to 
an entrance loss and fluid viscosity. 

Tables 13 through 16 contain the direct and cross-cou- 
pled stiffness (K, k) and damping coefficients (C, c) for 
the HPOTP seals. An examination of Table 13 indicates that 
the K values for seals 2 and 4 are significantly larger than 
for the other seals. This result is explained by longer 
lengths and larger AP's for these seals. For the remaining 
coefficients (k, C, c) , seal 2 is dominant, while the magni- 
tudes of coefficients for seal 4 are much reduced. The re- 
duction in magnitude of these coefficients for seal 4 is 
explained by the differences in the functions no and m and 
the dependence of the coefficients (k, C, c) on the transit 
time T defined by. 

T = £,/v = J.Y.A/W 
Y = specific weight 
A = flow area 
w = weight flow rate 

The ratio of y/w is larger for the lox of seals 1 and 2 than 
for the liquid hydrogen of seal 4. 


31 


CHAPTER III 

LINEAR STABILITY ANALYSIS FOR SSMF TURBOPUMPS 
3.1 INTRODUCTION 

The vibration or precessional motion of a rotor in typ- 
cially and ideally at its running speed, i.e., the speed of 
precession and speed of rotation of a rotor are equal. 

This mode of operation is customarily described as "synchro- 
nous". The amplitudes associated with synchronous motion 
are directly proportic t. J to the unbalance, and the custo- 
mary operational problem associated with synchronous motion 
is that of operation at or near a "critical speed". This 
latter condition arises when the running speed of a rotor 
coincides V7ith a rotor natural frequency, and the large am- 
plitudes associated v/ith critical speed operations are as- 
sociated with a resonance phenomenon. 

By contrast, many .rotors have experienced large-ampli- 
tude motion associated with unstable subs;. nchronous whirling. 
This type of motion typically begins at an operating speed 
above a critical speed. As the running speed is increased, 
the amplitude of subsynchronous motion increases, frequently 
leading to hardv/are failures. The characteristics of this 
unstable motion are independent of the imbalance, and can 
arise due to a variety of physical mechanisms. A partial 
listing of these mechanisms follows: 

(a) Hy.pteretic or internal friction damping due to shrink 
fits, friction in gear type couplings, or material 
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damping. 

(b) Hydrodynamic fluid film bearings and seals. 

(c) Aerodynamic cross coupling forces, 

(d) Dry friction whirl (rubbing) . 

An exr;ellent survey of rotor-bearing stability is provided 
in Ref. [12] by Gunter. 

Historically, flexible rotors supported on flexible 
bearings have been prone to unstable whirling modes. Speci- 
fically, rotors which are operated at speeds that are m.uch 
higher than rotor critical speeds are more likely to exper- 
ience unstable whirling motion. 3oth of the high pressure 
SSME turbopumps have top operating speeds which are consi- 
derably above their first critical speeds, and a>-e poten- 
tially subject to unstable whirling motion. This is parti- 
cularly true of the HPFTP, which employs a soft bearing 
support system. The analysis procedure developed in this 
study to examine the linear stability characteristics of the 
SSME turbopumps is discussed in the following section. 

3.2 LINEAR STABILITY Al^ALYSIS PROCEDURE 

The analysis procedure developed in this study is pro- 
bably similar to that outlined by Black [4]. However, a 
lack of detailed discussion in Black's paper m.akes a direct 
coiiiparison difficult. 

The linear stability analysis to be developed here will 
be based on the "conventional" modal Jeffcott model of Ap- 
pendix A, v;hich accounts for the structural dynamics of the 


33 


rotor-bearing system. In addition ^ the analysis v^ill ac- 
count for the seal-force definitions of the preceding chap' 
ter, speed-dependent properties of the ball bearings, and 
asymmetric bearing support stiffness. The appropriate re- 
sults from Appendix A are 
(9x> 

tPx) 

(Fj,) 

(Hjj) 

(^) 

The formulation of Eq. (11) directly accounts for gyrosco- 
pic cross-coupling via the [CM] matrix. The matrix [2S H 
is diagonal v?ith the elements and accounts for inter 

nal rotor damping. Since the stability of the rotor (as op 
posed to its response) is of interest, the modal imbalance 
vectors of Eg. (12) dropped. The modal 

force definition of Eq. (13) will be used to account for 
seal forces and. speed-dependent changes in the bearing char 
acteristics . Since no moments are involved, the vectors 
(My), (M^) can be dropped from Eq, (13). 

Incorporation of the effects of the seal-force defin- 
itions of the preceding chapter into the governing differ- 
ential Eq. (12) is simplified by defining the rotor's posi- 


= (Fy) - ^[CM](q^) - I2?ll{(gjj) + (^(qv)} 


“ [A] (q^j) + ^^(P^) + (P-y) 


X' 


(Fy) + <j>[CM](qjj) - [2eX]{(qy) 

- [AKg^) + 


« (q^) 3 


( 12 ) 


[Agl^(fx) + [Ag]'^(My) 

lAg]^(fy) - [Ag)^(My) 


(13) 


tAe]^(qx)r (Ry) = lAgl^(qy) 

= tAgl^(qy) 


(14) 


I 


tion and velocity at seal location as follows 
(«xs> = SSI (9x> ' = SS' 

(R^g) = [S] (4x) ^ 

The reduced eigenvector matrix [S] contains the rows of the 
[Ag] matrix which corresponds to seal locations, and is an 
nsxm matrix where ns is the number of seals and m is the 
number of modes retained for integration. 

Prom Eqs. (4) and (10) , the components of the seal 

force at the ith seal are defined by 

”'^XSi ^i^Si ^i^Si ^i^Si °i^Si 

"^YSi " '“^i^Si ^^i^Si " ^Asi *^i^Si 
and the vector of seal forces can be stated 

(f^g) = [K] (Ej^g) + [kl (Ryg) + [C] (R^g) + tc] t^g) 

(f^g) = **[kHR^g) + [k](Rys^ " ^ Ec]t^g) 

v;here [K] , [k] , [C] , and [c] are diagonal matrices with el- 

ements K, k, C, and c, respectively. The modal force defi- 
nition in terras of the seal forces is 

(F^) = tSl^(fjjs)/ (^Y^ == [Sl'^Cf^g) 

From Eqs. CIS) and (17), this result becomes 

_(p ) [SK] (q^) + tsk] (q„) + [SC] (q^j) + [Sc] (q^) 
XXI ^ ^ 

-(F ) =.-[Skl(qx) + [SK] (q^) - [Sc] (qj^) + [SC] (q^) 


where 

[SKd. = [S]^[K][S], [Sk] = [S]'^Ik][Sl 

[SC] = [S]^[C][S], [Sc] = [S]^[c][S] 

Substitution of Eq, (19) into Eq. (12) will yield the de- 
sired differential equations including the effects of seal 
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coefficients. 

A development similar to the above is used to account 
for speed-dependent changes in the bearing stiffnesses, and 
(v 7 here appropriate) changes in the bearing support structure, 
and yields a modal force definition of the form 

- (Fjj,) - tKXl (qjj) , - (F^) = tKY] (q^) (21) 

Substitution of Eqs. (19) and (21) into Eq. (12) yields the 
following homogeneous modal differential equatioiis 
(q^j) + {[2?X] + [SC]}(q^) +(HCM] + [Sc] Hq^) 

{[A] + [SK] + [KX]}(qx) + + [Sk] } (q^) = 0 

... (22) 
(g^) + {[2?X] + ESClHqy) " {<|>tCM].+ [Sc] } (q^j) 

+ {[A] + [SK] + [KY]}(g^) - {<^[2^X] + [Sk] Hq^) = 0 

These equations have the form 

(q) + [C] (q) -1- [K] (q) ^ 0 (23) 

where 

(q)^ = {(qx)"^/ (qy)^} 

The substitution (q) = (v) reduces the second order matrix 
differential Eg. (23) to the following first order differ- 
ential equation 



EK]! W 
0 


0 


(24) 


The eigenvalues of the coefficient matrix in Eq. (24) 
define the stability properties of the rotor. If m modes 
are retained in the analysis, the differential Bq. (24) is 
of 4m order and has 2m sets of complex conjugate eigenvalues; 
m of these sets correspond to forward processional motion. 
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and m sets correspond to backward processional motion. A 
complex-conjugate eigenvalue of the fom — ± indi- 

cates that the rotor is unstable, and that its motion will 
increase exponentially at the frequency 


3.3 AN APPROXIMATE LINEAR STABILITY ANALYSIS PROCEDURE 

Results from the above procedure have the advantage of 
accuracy in that the effects of all modes are accounted for; 
however, the results are difficult to interpret physically. 
Specifically, they indicate that the rotor is either stable 
or unstable, but do not directly provide insight into the 
physical mechanisms that either enhance or degrade stabil- 
ity. To gain a better appreciation for the physical mech- 
anisms , the following approximate stability analysis was 
employed , 

The assximption is made that the stability properties 
of the ith modal pair are independent of the re- 

maining modal coordinates. Hence, from Eq. (22), the gover- 
ning differential equations for this modal pair are 
g^i + {24 Xi + sc^^) q^j. + (|CM.^ + Sc^.)qy^ 

+ + SK^. + KX^i)gxi + ^2?. + Sk, .)qy. = 0 

C25) 

. U.2 + + KY. = 0 

The assumption that the modal pair q„. r Ov are uncoupled 

AJL xX 


from the remaining modal coordinates has been enforced by 
retaining only the diagonal elements of the matrices [SK] , 


[Sk] , [SC], [Sc], [BK] , and [CM]. 


The characteristic equation for the differential Eq. 
(25) has the form 

s'* -t- a is^ + ^ 2 ®^ ^ 3 ^ + at,. = 0 (26) 

The Routh-Hurwitz conditions for system stability require 
that the coefficients of this equation a^^ be greater than 
zero, and that the following tv7o conditions be satisfied 


A 2 - 1 

a 1 1 

O 

A 

A 3 ” 

ai 

1 

0 


a.3 a£ 



^3 

a£ 






0 

ai. 

as 


The coefficients a^ are always positive, and the require- 
ments from condition { 21 ) are 

,(X| H- SK^. + K.) + (25. + SC^.)2 

(2B) 

+ (0CM. . + SC^. ) (|CM. . + SC. . ~ > 0 


+ (X? + K, + SK . ) + H.UCM.. + SC.. - 


XX 


H. ) > 0 


(29) 


where 


2K. = KX . . + KY . . , 26 K . = KX . . - KY. . 

X XX XX X XX XX 


(30) 


H. = (2c^X.J + Sk^.)/(25iX,^ + SCj^) 

The stability requirement of (28) and (29) are always sat- 
isfied if 


^CM. . + Sc. . - H. > 0 

^ XX XX X 


(31) 


However, for the SSME turbopumps, this condition is normal- 
ly violated. This point is better appreciated if one notes 
from Eq. (20) that the modal seal coefficients Sk.., SC.. 



are proportional, respectively, to the physical seal coef- 


ficients k., c.. Prom Eqs. (4), (5), and (ID) the physical 

X X 

seal coefficients k^, c^ are in turn proportional, respec- 
tively, to vi^T/2 and Hence, the coefficient Sk. . is 

proportional to while is not, and the definition of 

provided by (29) implies that is also proportional to 
i . Condition (31) can then be restated as 

Sc.. + i (CM. . - h ) H. = ih. (32) 

For the SSME turbopumps, hj|^ > which implies that the 

favorable gyroscopic stiffening effects for these rotors is 
less pronounced than the combined detrimental effects of in- 
ternal rotor damping and the seal cross-coupling stiffness 
coefficients . 

Given that condition (32) is violated as ^ is increased 
condition (29) will normally be violated prior to condition 
(28) . Hence subsequent considerations V7ill be restricted 


to condition (29) , which can be restated as 

(6K.)2/{ (2c. a. + sc. .)2 + (^CM. . + Sc . . ) ^ } 

' x' ' X X XX ^ XX XX (33) 

+ ( A? + SK. . + K. ) + |h.CM, , - (h. - CM. . ) > 0 

The term |h^Sc^j^ is insignificant in this relationship be- 
cause (a) the modal seal *'r*'A: tfficient Sc^^^ is comparatively 
small (i.e., SK. . > Sk. . > SC.. > Sc..), and (b) this term 

XX XX XX XX 

is proportional to | , while the right-hand term is propor- 
tional to The first. ;|^rm in (33) arises due to asymme- 

try in bearing support s; J.ffness and vanishes for a "typi- 
cal" symmetrically supported rotor. The term (A? + SK^^^+K^) 

page is 
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represents the axisymmetric stiffness of the rotor-bearing- 
seal system. The larger this term, the larger ^ must be to 
cause an onset of instability, which explains in part why 
stiffer rotors tend to be more stable than flexible rotors, 
and why instabilities are normally associated with the 
first or lowest-natural-frequency mode. 

The stability criteria (33) includes the effect of gy- 
roscopic stiffening, vrhich was neglected by Gunter and 
Trumpler [13] in a previous analysis of the effect of stiff- 
ness asymmetry, lor large condition (33) reduces to 

,)2 + ! ^ + K^)- }2h. (h^ - CM^.) > 0 

which indicates that the gyroscopic coupling term reduces 

the favorable effect of stiffness asymmetry as running speed 
I is increased. 


3.4 LINEAli STABILITY ANALYSIS FOR THE HPFTP 

The analysis procedures of the preceding section 
were employed to establish the stability characteristics 
of the SSME turbopumps. The rotordynamic models used for 
the turbopumps are discussed in References [1] and [2] , 
with sujjplementary data and discussion provided in Appen- 
dix B. 

As noted in the preceding section, the t 3 igenvalues of 
the coefficient matrix in Eq. (24) define the linear stab- 
ility characteristics of the rotor model. The speed-de- 
pendent terms in this matrix arise from the seal coeffi*- 
cient matrices [SK] , [Sk] , [SC] , [Sc] , and the speed-depen- 
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dent changes in the bearing stiffnesses. Changes in the 
bearings stiffness as a function of speed were obtained di- 
rectly from Eg. (B.2) , and the seal coefficients of Tables 
5 through 8 were least-square curvefitted by cubic polyno- 
mials in I for the speed-set of Eg, (9) . Specifically, cu- 
bic polynomial coefficients were calculated for the coeffi- 
cients Cj^ for each seal- A cubic polynomial 

did a satisfactory job of fitting the coefficients of most 
seals , particularly the dominant seals 3 and 6 . 

A computer program was developed to assemble and eval- 
uate the coefficient matrix of Eg. (24 ) , and to subsequent- 
ly calculate its eigenvalues. The eigenvalue calculation 
was performed via a QR algorithm using the I.B.M. Scienti- 
fic Subroutine Package subroutines HSBG and ATEIG, 

The results developed in this section demonstrate the 
dependence of the HPFTP rotor stability on variations in 
(a) the bearing support stiffness and (b) internal rotor 
damping. As noted in Appendix B, the HPFTP bearing support 
has a ’’design" radial stiffness of 

kgp(pump) = kg^ (turbine) = kg = 2,17 x 10^ Ib/in, (34) 

which yields critical speeds at approximately 9,000, 17,500, 
and 48,000 rpm. The calculation of these critical speeds 
includes the speed-dependent properties of the bearings and 
seals . 

Tests of the bearing support structure hardware have 
yielded the following load-deflection data. 
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load (lb) 


deflection (in) 

2.6 X 10"3 
3.75 X 10"3 
4.5 X 10"3 
5.0 X 10“^ 


(Ib/in) 

1.38 X IQS 
1,74 X 10® 
2.67 X 10® 
4.0 X 10® 


The "stiffness" kg is defined as the change in load divi- 
ded by the change in deflection, and is the local slope of 
the load-deflection curve. Deflections greater than .005 
in. would exceed available clearances, and are not of in- 
terest. The data indicates a significant nonlinear har- 
dening of the bearing support structure as deflection mag- 
nitudes are increased. However, this effect is not pro- 
nounced over the expected range of deflection magnitudes, 
which for stable operations would normally be less than 
-003 in. From the above data and considerations, the 
bearing support stiffness 

k = 1.4 X 10® Ib/in (35) 

was selected as more representative of the actual hardware, 
and a reasonable lower-bound choice. The calculated criti- 
cal speeds for this support stiffness are 8,000, 15,200 and 

■r 

45,500 rpm. 

The above cited radial bearing stiffnesses are all 
ax i symmetric following the conventional practice in turbo- 
pump design. However, the stability criterion (33) demon- 
strates that rotor stability can be enhanced by introducing 
asymmetry into the bearing support structure. This ap- 
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proach has previously been analyzed by Gunter and Trump- 
ler [13 ] , and was investigated here for the following asym- 
metric support stiffness selection 


kgpjj - 2.3 K 10® Ib/in ^stX ” 2.3 x 10® Ib/in 
kgpY = 1.38 X 10® Ib/in ~ 1.38 x 10® Ib/in 


(36) 


The calculated critical speeds for this set of bearing 
stiffnesses are approximately 9,100, 17,300, and 50,000 rpm. 
The bearing stiffness choice in Eg. (36) was obtained by 
holding ^gpyf ^sty constant, and increasing ^stX 

til the second critical reached its cited value. An addi- 
tional increase in kg^^^ is not desirable, since it 

would place the second critical speed too near the operating 


range . 

The analysis of the preceding section demonstrates 
that internal rotor damping, as defined by the modal dam- 
ping factors acts to destabilize a rotor as (|> is in- 
creased. Hence the results of stability analysis are strong' 
ly dependent on the choice for and for the HPFTP are de- 
pendent on the first two modal damping factors 5 i, ?2- From 
past test and operational experience. Rocket dyne personnel 
suggest upper bounds for ?i, ?2 order of .01 to .02; 

hoviever, (as always) the correct value for these factors 
are highly urcertain. In this section, the choices for in- 
ternal damping to be considered are (a) = ?2 “ Of 

~ ^2 “ *01/ with ?3 = = .02. The true situation is 

assumed to be bracketed by these values. 


The results of stability analysis for the HPPTP rotor 
with ^ 1 = ?2 = 0 axisymmetric stiffnesses of Eqs. 

(34) and (35) are presented in figure 4(a), and consist of 
the roots associated with the first tv/o modes at the run- 
ning speed set of Eq- (9) . The results are presented in 
root-locus format in figure 4 (a) . Note that the roots are 
the negative of the eigenvalues? hence, stable and unstable 
roots are located in the left and right half planes, respec- 
tively. The upper and lower half planes are mirror images? 
hence, only the upper half plane is illustrated. The data 
marks on the curves denote the speeds of Eq. (9) . The re- 
sults indicate that the second mode is stable over the oper- 
ation sj-.eed range, but the onset speeds of instability for 
the first mode at kg = 2.17 x 10^ Ib/in and kg = 1.4 x 10® 
Ib/in are 19, 500. and 16,800 rpm, respectively. The predic- 
ted whirl frequencies (the frequency at which the rotor 
would whirl subsynchronously) are 9,850 rpm (kg = 2.17 x 
10® Ib/in) and 8,500 rpm (kg = 1.4 x 10® Ib/in). 

Figure 4 (b) illustrates the same results as 4 (a) ex- 
cept that Cl = Cz = .01, In this case the first mode is un- 
stable over the entire speed range of Eq. (9), and the onset 
speeds of instability for the second mode are 22,900 rpm 
(kg = 2.17 X 10® Ib/in) and 19,200 rpm (kg = 1.4 x 10® Ib/in) 
An extrapolation of the results of figure 4 (b) yields pre- 
dicted onset speeds of instability for the first mode at 
12,200 rpm (k„ = 2.17 x 10® Ib/in) and 11,000 rpm (k- = 


1.4 



J 



. figure 4a 

EEETP st^ility analysis results for = C 2 = 0*0 and the speed range 

j- C 

of Bq* (9) with (a) ~ 2.17 x 10 Ih/in -j4 — / and (b) 1«4 x 10 Ib/in 



•/O to ZO 30 40 


figure 4b 

HE3BTP Stability analysis results for ? 1 *= C a “ 0* 01 and the speed range 

of Eq. (9) with (a) k = 2.17 x 10^ Ib/in ~ir£r- , and (b) 1.4 x 10^ Ib/in 
# s 
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The results of figures 4 (a) and 4 (b) v/ere checked by 
using the approximate stability criterion (33) of Section 
3.4 for the first mode. The approximate procedure agreed 
in all cases with the more exact procedures of Section 3.3. 

Note was made in the preceding chapter that the stiff- 
ness and damping coefficients for seals 3 and 6 are an or- 
der of magnitude larger than the coefficients for the other 
seals. The above stability analyses were accordingly repea- 
ted with the coefficients of all seals except seals 3 and 6 
set to zero, with no appreciable change in the results. 

The results of figures 4 (a) and 4 (b) support the fol- 
lowing general conclusions: 

(a) The HPFTP rotor is predicted to be unstable because of 
the seals alone (no internal damping) . 

(b) Rotor stability is extremely sensitive to internal mo- 
dal damping factors ? i, tz, 

(c) Feasible increases in the axisyrametric radial stiffness 
k„ does not have an appreciable influence on rotor sta- 
bility . 

Figure 4 (c) illustrates the results of stability analy- 
sis for the asymmetric support stiffness set of Eg. (36) 
with = ?2 == 0^ ^nd Cl = 5?. = .01. For = 0^ the 

second mode is stable for the speed range of Eg. (9) , and 
the onset speed of instability for the first mode is 32,600 
rpm. For c i = Ca " *01r the onset speeds of instability for 
the first and second modes are, respectively, 30,300 rpm and 
26,300 rpm. Hence, one has the surprising result that in- 



figure 4c 

HEPTP stability analysis resu3.ts for the assOTnetric support stiffness 
choice of Eq. (35) and the speed range of Eq. (9) vnLth (a) = Cz = ' 


* 


and (b) 51 = 52 = *01 — 


t 
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ternal damping, while reducing both onset speeds of instab- 
ility, yields a lower onset speed of instability for the 
second mode than the first. The explanation for this re- 
sult is that CM 11 is an order of magnitude smaller than 
CM 22 » hence, from condition (33), the stabilising effect of 
asymmetry is more pronounced for the first mode than for 
the second. A comparison of figure 4(c) with figures 4(a) 
cind 4 (b) demonstrates that the introduction of support asym- 
metry can markedly improve the stability of the HPPTP rotor. 
However even with the maximum support asymmetry provided by 
Eg. (36) , a rotor instability continues to be predicted. 

The nonlinearity of the bearing support structure has 
the potential of providing some "natural" asymmetry with- 
out modification of the bearing carriers. Specifically, a 
steady-state deflection of the rotor in the X-Z plane due 
to the hydrodynamic side load causes the support structure 
to be stiff er in the X-Z plane than in the Y-Z plane. Since 
the steady state displacement of the turbine bearing is sig- 
nificantly larger than the pump bearing, this nonlinear har- 
dening effect would be felt primarily at the turbine bearing. 
From these considerations, stability analysis was performed 
for the HPFTP rotor for the asymmetric support stiffness set 


SpX “ 

X 10^ 

Ib/in 

■'stx = 

X 10^ 

Ib/in 

SPY = 

X 105 

Ib/in 

kg„ = 1.38 

X 10^ 

Ib/in 


This degree of asymmetry had no appreciable influence on ro- 
tor stability. 


The results of figures 4(a), (b) are not encouraging, 
since they generally predict a sub synchronous rotor insta- 
bility, which has the possibility of serious damage to the 
turbopump. From past experience, the bleak situation of 
figure 4(b) seems very unlikely; however, the situation 
illustrated in figure 4(a) is consistent with past rotor 
instability experience [15] . Specifically, the predicted 
onset speeds of instability are approximately twice their 
associated whirl speeds. From a percentage of critical 
damping viewpoint, the rotor is never "severely" unstable, 
e.g., at FPL the percentage of (negative) critical damping 
in figure 4(a) is on the order of 0.5. The smallness of 
this value will be immaterial if tests prove the rotor to 
be unstable; however, it does suggest that small changes in 
the rotor dynamic model could yield a contrary prediction; 
viz . , a stable instead of an unstable rotor . 

The principal point of uncertainty in the model which 
yields a prediction of instability is the stiffness and dam- 
ping coefficients for the seals. The seal models developed 
by Black are for a centered position of the shaft relative 
to the seal. For the HPFTP rotor, the combination of a 
large side load and soft bearing mounts insures violation 
of this condition. Black has performed nonlinear calcula- 
tions (but not tests) which suggest that his linear seal mo- 
del is reasonable for eccentricity ratios on the order of 
0,5; however, the validity of Black’s nonlinear "turbulent 
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bulk flov?" model can also be questioned* The justifica- 
tions for using Black's model are (a) it appears to be rea- 
sonable and for a centered position is supported by test 
data, and (b) it is the only model presently available. 

In addition to the seal stiffness and damping coeffi- 
cients and internal rotoi~ damping, rotor stability probably 
depends on the follov/ing additional physical phenomena: 

(a) Ball bearing damping. The fact that loaded ball bear- 
ings operating at high speeds generate a great deal of 
heat suggests that ball bearings would be a source of 
external rotor damping. However, tests of rotors on 
ball bearings have indicated little or no damping, and 
a zero-damping assumption is generally followed in ro- 
tordynamic analysis 

(b) Small radial clearances are provided at the bearings on 
both turbopumps to allow axial motion of the rotor. 
These clearances are obviously filled with the working 
fluid of the pumps, and the annular surfaces associated 
with these clearances could act like squeeze-film dam- 
pers , which would in general improve the stability pro- 
perties of the rotors. No consideration was made of 
these clearances in the present study. 

(c) Aerodynamic Cross Coupling. Alford [14] has proposed 
an aerodynamic cross-coupling mechanism as the cause of 
some rotor instabilities. The physical explanation for 
this instability is the difference in efficiency of 
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turbine blades in the vicinity of r-educed or enlarged 
clearances due to rotor displacements. However r very 
little progress has been made since Alford’s original 
paper in apiiori calculation of aerodynamic cross-cou- 
pling coefficients, and as things now stand this mech- 
anism tends to be invoked to explain rotor instabili- 
ties which depend on power levels as well as 
speed. Presumably, the turbine wheels for the SSME 
turbopuraps could introduce some aerodynamic cross-cou- 
pling which would degrade the rotor's stability char- 
acteristics; however, no consideration was made in this 
study of aerodynamic cross-coupling. (Note is made, 
however, that the analysis procedures of Sections 3.2 
and 3 . 3 can readily account for aerodynamic cross-cou- 
pling coefficients.) 

If the rotor instability- predicted by figures 4a, b is 
encountered during the SSME test program, the easiest and 
most direct corrective modification to the turbopump would 
be the introduction of stiffness asymmetry as illustrated 
in figure 4c. However, a judgement on the probable success 
of this approach would have to be based on the measured 
characteristics of a rotor instability, e.g., onset speed 
of instability, v/hirl frequency, rate of divergence etc- 
In the following chapter, an alteriiative (or supplementary) 
approach for the improvement of rotor stability is consi- 
dered, which involves the attempted development of squeeze- 
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film damper forces at the annular surface between the bear- 
ing cari*iers and the housing of the turbopiimps. 

The opinion of this investigator is that rotor stabil- 
ity is not likely to be significantly improved by external 
coulomb damping devices or by feasible modifications of the 
seals in the HPFTP. Ccalomb damping is not generally effec- 
tive after an instability has been initiated, and seals 3 
and 6 presently have configurations 'Thich are close to those 
recommended by Black [4] . Rotor stability could possibly 
be improved by a program to reduce internal rotor damping 
by minimizing rubbing. This approach has been followed in 
dealing v>ith shrink-fitted wheels on steam turbines [15] . 

3.5 LINEAR STABILITY ANALYSIS RESULTS FOR THE HPOTP 

The analysis performed in the preceding section was 
basically followed here in analyzing the KPOTP rotordyna- 
mic stability. The HPOTP rotordynamic model is reviewed 
in Appendix B, The requisite seal coefficients are pro- 
vided in Tables 13 through 16, and were curve-fitted via 
cubic polynomials in | . The definition of speed-depen- 
dent changes in the bearing stiffnesses were obtained di- 
rectly from Eqs, (E.4) and {B.5). 

From a rotordynamics viewpoint, the principal dif- 
ference between the HPOTP and the HPFTP is that the HPOTP 
does not employ a flexible bearing carrier. With the ex- 
ception of small radial clearances, vfhich are provided to 
permit axial motion, the bearings are supported directly 


by the turbopiimp housing. By the nature of this simpli- 
fied design, tliere is little that can be done to this tur- 
bopump to either enhance or degrade its rotordynamic sta- 
bility. 

The first tv;o critical speeds of the HPOTP have been 
calculated by Rocketdyne to lie at approximately 13,000 
and 40,000 rpm. Hence, only the first critical speed lies 
v/ithin the operating range, and only the first forx^rard 
precessional mode is potentially unstable. The next three 
modes are comfortably stable for all conditions examined 
here, and the balance of this discussion will deal exclu- 
sively with the first mode. 

Figure 5 (a) illustrates the results for the first 
mode of stability analysis v/ith = 0.0 ^ver the speed 
range of Eg, (11). The results indicate that the rotor 
is stable over the speeds considered when acted on solely 
by the seal forces. The predicted first critical speed 
associated with the results of figure 5(a) is 14,600 rpm. 
An explanation for the fact that this value is higher than 
the 13,000 rpm prediction of Rocketdyne is provided by 
the comparatively large values of the direct stiffness co- 
efficients in Table 13. The stiffness at seal 4 is of 
particular importance in this regard, since the first- 
mode eigenvector (Table B.2, Appendix B) has comparatively 
large deflections associated with the over-hung turbine 
wheels at this seal location. Hence, the comparatively 
large direct stiffness at seal 4 location elevates the 
critical speed. 


1 
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Figure 5 (b) illustrates the results for the first 
mode of stability analysis with = .02 over the speed 
range of Eq. (11) . The results predict an onset speed of 
instability at approximately 15,000 rpm, again with the 
first critical speed at 14,600 rpm. The results of this 
figure are probably unrealistic, since rotor instabilities 
are more frequently associated with speeds that are on 
the order of twi< ' a rotor natural frequency, whereas 
this predicted instability speed is only slightly greater 
than its associated critical speed. 

The results of figures 5(a) and (b) were checked and 
confirmed by the stability criterion (33) of Section 3.3. 

A comparison of figures 5 (a) and (b) emphasizes the 
dependence of stability on the internal modal damping 
factor ^ 1 . An inspection of Table B.l of Appendix B de- 
monstrates that rhe HPOTP first mode is a bending mode as 
opposed to the first two HPFTP modes which have a rigid- 
body character. Hence, a comparatively large modal dam- 
ping factor might be anticipated. However, the basic 
structural element in the HPOTP rotor is a solid shaft, 
v^/hile the HPFTP rotor is basically a group of wheels 
held together by a throughbolt. The parts v/hich are con- 
centrically fitted to the outside of the HPOTP rotor out- 
board of the turbine bearings ai"e not designed to carry 
a bending load, and the rubbing associated with their in- 
terference fit is presumably small. These factors argue 
for a comparatively small t,\, with a consequent improve- 
ment in the stability picture. 
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figure 5b 

HEOTP stability analysis for ?i = 0.02 and the speed range of Bq. (11) 



figure 6 

Recjcrmiended squeeze-film daitper design from reference (16) . 

PAGE IS. 

Of fODE QUAIM’ 



55 


The linear stability analysis results of this and the 
preceding section were generated on an I.B.M. 360-165 compu- 
ter. The computer time requirements for the speed set of 
Eg. (11) used in analyzing the HPOTP were approximately 8 
or 9 seconds, with a comparable time requirement for the 
HPFTP speed set of Eq- (9) , 

3.6 SUMMARY 

a review of the stability results for the HPFTP and 
HPOTP supports the following general conclusions. 

(a) The HPFTP is probably going to encounter stability 
problems within its range of operations; however, if 
an instability is encountered, various procedures can 
be attempted to improve its stability properties. 

(b) The HPOTP is inherently more stable than the HPFTP, 
and is less likely to experience stability problems 
over its range of operations. However, if this unit 
does have an instability problem, there are a limited 
number of obvious approaches for improving its stab- 
ility characteristics. 
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CHAPTER IV 

TRANSIENT SIMULATION OF SSME ROTORDYNAMICS 
4,1 INTRODUCTION 

As noted in the introduction, prior studies by this 
investigator of SSME rotor dynamics [2,3] employed tran~ 
sient simulations to predict amplitudes of rotor motion 
during critical-speed transitions and at FPL. In these 
investigations, sufficient external damping was provided 
to insure stability. The results of the preceding sec- 
tion would suggest that this was an overly optimistic 
approach, since both turbopumps are predicted to be either 
veiry lightly damped or unstable at full-power levels . If 
the results of Chapter III are in error, and the rotors 
are in fact stable at full-pov/er levels, the predicted 
amplitudes of Ref ^ [2-3] should be reasonable, since the 
PPL speeds are well removed from critical-speed locations, 
and are comparatively insensitive to rotor damping. How*^ 
ever, if the rotor damping is as low as predicted by the 
results of Chapter III, the amplitudes of the rotors 
during critical-speed transitions are going to be signi- 
ficantly larger than predicted in Ref. [2,3]. In particu- 
lar, the predicted rubbing problem for the HPOTP turbine 
floating-ring seal during shut-down critical-speed tran- 
sition v/ill be more severe than previously suggested. 

The following results of Chapter III, with regard to 
rotor stability, established the objectives of transient 
rotordynamic simulations for the present study. 
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(a) Very serious stability problems are predicted for the 
HPFTP, and several alternatives can be considered to 
improve the situation. 

(b) Stability problems are possible for the HPOTP^ but 
less likely than for the HPFTP. Moreover, there are 
few obvious options presently available to improve 
the stability characteristics of the HPOTP. 

In view of these results, transient rotordynamic in- 
vestigations of this study were concentrated 
primarily on the application of squeeze-f ixm dampers to 
improving the stability characteristics of the HPFTP. A 
proposed squeeze-film damper design for the HPFTP is dis- 
cussed in the following section. 

4.2 SQUEEZE FILM DAMPER DESIGNS FOR THE HPFTP 

Figure 6 is taken from Ref. {16] and illustrates a 
"typical” squeeze-f i].m damper installation. In this de- 
sign, the outer race of a ball bearing is attached to a 
flexibly supported non-rotating cylinder. A squeeze-film 
damper is contained in the annulus at the outer surface 
of this cylinder. The fluid is a comparatively viscous 
oil, which is both pumped into and confined v/ithin the 
annulus. The damping forces developed via a squeeze f'? ■ 

damper differ from those developed at a journal bearing 
because the ball-bearing housing/damper journal is allov/ed 
to move radially in response to the rotor's motion, but 


does not rotate. 
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The design concept illustrated in figure 6 is clear- 
ly not compatible v;ith a liquid oxygen or liquid hydro- 
gen turbopump* Specifically, the problems involved in 
supplying a separate oil supply in a cryogenic environ- 
ment would be staggering. However, the design of the 
bearing carrier structure in the HPFTP is such that its 
outer surface executes tjie desired motion of a squeeze 


film damper, i.e,, i i '.;., esses in response to the bearing 

reaction loads , but does not rotate , Hence , a modifica- 
tion of the annulus at the outer surface of the bearing 
carrier structure as illustrated in figure 7 could yield 
squeeze film damper surfaces at both bearings in the HPFTP. 
These designs are certainly not as attractive as the 
ideal design of figure 6 , since the purap-and turbine-end 
bearings will contain liquid and gaseous hydrogen, as 
opposed to a comparatively high viscosity oil. Moreover, 
the hydrogen in these dampers is not pres sure- supplied 
through a central groove, nor is its axial flow con- 
strained by piston rings. These design modifications 
have the principal advantage of being compatible with the 
present turbopump desigi?'. In the balance of this chap- 
ter, the assumption is made that the proposed design mo- 
difications of figure 7 have an adequate supply of working 
fluid. 

The modifications figure 7 use a constant-clear- 
ance "conventional" design. An unconventional "partial" 
squeeze-film damper design v;as also examined in this stu- 
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figure 7 

Proposed squeeze-fi3m damper design nodifications for the HEETP 
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dy. This design utilizes a much-enlarged radial clear- 
ance over half of its circumference. The stability re- 
sults were indistinguishable for the conventional and par- 
tial squeeze-film-damper designs; however, an adequate 
supply of working fluid would be more easily assured with 
the partial design. 


4.3 HPFTP SQUEEZE FILM DAMPER MODELING 

The results of (the excellent) reference [17] by 
Kirk and Gunter was used here to model the nonlinear load- 
deflection characteristics of a squeeze-film damper. 

Kirk and Gunter use the short-bearing approximat ^on to 
obtain the follov/ing definition of the force components 
acting on a rotor from a squeeze film damper of length £ , 
radial clearance 6 , radius r, and fluid viscosity p 




S0 


d(£r0) 


(37) 


In this definition for the force components, the average 
pressure at angular position , defined by 


P(8) = 

(fi - Rjj.C6 - R^S0)3 


(38) 


is integrated over the circumference of the squeeze film. 
In Eq. (38), c9 = cosO, s9 = sin6 , and the position and 
velocity of the squeeze-film damper cylinder are defined 
respectively by Ry and Rjj, Ry. The term squeeze-film 

damper cylinder refers 
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here to the ball-bearing housing/damper jouj.'nal of figure 
6 or the ball-bearing carrier of figure 7. The effect of 
cavitation is enforced in reference 117] by requiring that 
the average pressure P{6) equal or exceed zero. Xn the 
present situation the pump-end damper contains gaseous hy- 
drogen at very high pressure, and the turbine-end damper 
contains liquid hydrogen at approximately 160 psi. (FPL 
condition) . Hence, the assumption is made that neither dam 
per will cavitate. Kirk and Gunter's recommendation of a 
Newton-Cotes quadrature formula for the evaluation of the 
integral in Eq. (37) was followed here, with an integration 
step of (tt/ 30 ) radians. 

Prom a transient simulation viewpoint, the squeeze- 
film damper reaction definition of Eq, (37) does not ba- 
sically introduce any new reaction-definition requirements 
beyond those normally required for a ball-bearing. Spe- 
cifically, the position and velocity of the damper cylin- 
der is required in Eq. (37) , and this is also the normal 
requirement for a bearing reaction definition. The pump- 
and turbine-bearing carrier eigenvector entries are, re- 
spectively, rows 28 and 29 of Table B.l, Appendix B, As 
noted in Appendix A, the correct definition of the posi- 
tion of the bearing carrier structures required the use 
of a residual correction procedure in this study. 

The squeeze-film damper reaction definitions of Eq. 

(37) can be used to obtain direct and cross-coupled stiff- 
ness and damping coefficients similar to those of Eq. (4) 
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for feals. Alternatively, for the purposes of design, "equi- 
valent" direct stiffness and damping coefficients can be cal- 
culated as a function of eccentricity ratio e. The physical 
and dimensionless stiffness (K^, K^) and damping (R^, R^) 
coefficients are defined in terras of the eccentricity ratio 
E by 


= 2e/(1-e2)2 

R^ = = ir/2 (1 -e2) 


(39) 


The dimensionless coefficients are obviously nonlinear 
functions of the eccentricity ratio , which approach an 
unbounded limit as the eccentricity ratio e approaches 
unity. 

The unsuitability of hydrogen as a squeeze-film 
damper fluid can be appreciated by calculating these co- 
efficients for an 'eccentricity ratio of c = 0,4, which, 
yields from Eq. (39), = 1,13? R^ = 2,04, The remai-^ 

ning data for the HPFTP at FPL are: 

I = 3930 rad/sec, r = 1,843 in, 

yCpump end) =* 2,96 x 10”2ibf sec/ft^ = 2.06 x lO^^lbf sec/in^ 

y (turbine end) = 2.18 x 10**2ibf sec/ft^ = 

1.52 X 10"“^lbf sec/in2 

£ (pump end) = £ (turbine end) = 1 in 
Prom these data, the follov^ing physical stiffness and 
damping coefficients can be calculated 


6 (in) 

R^db sec/in) 

K^(lb/in) 

.001 

4.21 

16,600. 

.002 

,526 

2,070. 

.003 

.150 

615 . 


These values are small in comparison to the seal stiffness 
and damping coefficients of Tables 5 and 1 , but from Eq. 
(39) , they increase sharply as the eccentricity ratio in- 
creases. 

4.4 TRANSIENT SIMULATION OF HPFTP MOTION WITH SQUEEZE-FILM 
DAMPERS AND AXISYmETRIC STIFFNESS SUPPORTS 

The results of Chapter III demonstrate that the HPFTP 
rotor becomes progressively less stable as the running 
speed I is increased; hence, all transient simulations for 
the purposes of examining rotor stability were conducted 
at FPL running speed; | = 37,470 rpm = 3930 rad/sec. The 
data of Eq. (40) was used for all simulations. The propo- 
sed damper redesign of figure 7 includes a modification of 
the ball-bearing carrier to achieve the damper length of 
one inch given in Eg. (40) . Maximization of the damper 
length is obviously desirable, since the damper forces 
are proportional to the length cubed. 

The transient simulation model was first employed to 
verify the stability results of Chapter IIIc The results 
of the transient simulation model agreed in all respects 
v;ith the linear stability analysis results of Chapter III, 
The simulation model was then employed to determine the 
effectiveness of the squeeze-film damper designs of figure 
7 in stabilizing the rotor. The capabilities of the si- 
mulation model were required to account for the signifi- 
cant nonlinearities of the squeeze-film damper force de- 
finition in Eq, (37) . 
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The imbalance used in the simulation model is as fol- 
lows 

1st stage impeller, 2.86 gm in 
2nd stage impeller, 2.96 gra in 

3rd stage impeller, 3.52 gm in (41) 

2nd stage turbine > 2.74 gm in 

1st stage turbine, 2.88 gm in 

The eccentricity distribution used to yield these imbal- 
ances was 

1st stage impeller, a^^ = 3.57 x lO~*^in 

2nd stage impeller, = 3.68 x lO^^in 

3rd stage impeller, a^ = 4,02 x l0”*^in (42) 

2nd stage turbine, a^ = 3.05 x 10~‘*in 

1st stage turbine, = -3.03 x 10~**in 

The data of Eq. (41) was provided by Mr. B. Rov/an, Rocket- 
dyne, and includes tolerance stackup and reassembly im- 
balance, The imbalance distribution of Eq. (42) was se- 
lected to be representative, as opposed to a worst-case si- 
tuation . 

In addition to the imbalance distribution of Eq. (42) , 
the HPFTP rotor is also subject to a hydrodynam'ic side 
load at the third-stage impeller discharge. The direc- 
tion of this load is fixed with respect to the pump hou- 
sing, and the load is defined here by 

fj, = -6.721 X 10~5|2 (43) 

v?here | is in rad/sec. This fixed direction side load is 
large in comparison to the imbalance loading of Eq. (42) . 

From a design vie^irpoint, the only variable which can 
be modified in tlie squeeze-film damper design is the clear- 
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ance. The results of Eq. (39j in t>e preceding section in- 
dicate that the clearance of a squeeze-film damper 
should be minimized. The procedure follov/ed in deter- 
mining feasible minima for the pump- and turbine-end 
dampers consisted of running the simulation model with 
sufficient external damping to ensure stable synchronous 
motion and then simply picking clearances that are suf- 
ficient to enclose the stable rotor orbits at the bear- 
ings. 

Applying the above procedure to the situation of 
figure 4(a) (kgp = = 2.17 x 10^ Ib/in; iJi = ^2 = 0.0) 

yields the following squeeze-film damper clearances 
6 (pump) = .002 in 

P C44) 

(turbine) == ,003 in 

The frames of figure 8 indicate the results of a simula- 
tion run with these clearances for the squeeze-film dam- 
pers of figure 7 . The initial conditions for this run are 
the final conditions of a run which was stabilized by ad- 
ding external damping; however, the only damping acting 
on the rotor in figure 8 is that provided by the seals 
and the two squeeze-film dampers- 

Figures 8 (a) and (b) illustrate respectively the ro- 
tor-fixed modal coordinates q^2 ’ Note that q^^^ is 

approximately an order of magnitude larger than qjj-2* 
high-frequency oscillations in these modal coordinates 
are caused by the hydrodynamic side load of Eq. (43) , If 
the rotor were excited solely by imbalance these variables 





figure 8e 

Eotor motion at the turbine-end squeeze-filra daipper in tlie X-Z plane 
for k„ = Z.17 K 10^ Ib/in; = Ca = 0. 
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figure 8f 

Botor orbit at the turbine^end. squeeze~ f i 3 . n i danger for — 2.17 x 10 Ib/in 

Cx = ?2 = 0. 



figure 8g 

Pun^i-end squeeze-film daiiper cylinder motion iix the X-Z plane for 
k_ = 2.17 X 10^ Ib/on; Ci = 52 = 0. 
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figure 8h 

Puinp“end. sgueeze~film dairper sylinder orbit for = 2.17 x iO Ib/in, 
Cl = C2 = 0. 

c. 



figure 8i 

■jiarbiue-end sq*j.eeze-film dtimper cylinder notion in the X-Z plane for 
k = 2.17 X 10^ Ib/in; Ci = Cz = 0. 
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figure 8j 

5 

Turbine^! jd sgaeeze-film danger cylinder orbit for k = 2.17 x 10 Ib/in; 
?1 = ?2 = 0 * 



figure 8k 

Pump-end squeeze-film dairgsar reaction n^agn.itude for ]c 


= 2.17 X 10^ 


Ib/in 
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would be constant during synchronous precession of rotor. 
The results of figure 8 (a) indicates that the squeeze- 
film dampers v/ith the clearances of Eq, (42) have complete- 
ly stabilized the first mode. However, figure 8(b) illus- 
trates that the second mode has a sub synchronous oscilla- 
tion which has limit cycled. Without the squeeze-film 
dampers, the first mode is unstable, V7hile the second mode 
is stable. 

Figures 8 (c) and (d) illustrate the motion of the ro- 
tor at the pump bearing. Figure 8(c) demonstrates the mo- 
tion of the rotor in the stationary X-Z plane, and 8(d) 
illustrates the rotor orbital motion. The static dis- 
placement of the rotor illustrated in these figures is 
caused by the hydrodynamic side load of Eq. (43), These 
figures illustrate that the pump ^nd of the rotor is pre- 
cessing synchronously. 

Figures 8(e) and (f) are similar to figures 8(c) and 
(d) , except that they illustrate rotor motion at the tur- 
bine-end bearing. Both of these figures indicate the pre- 
sence of a subsynchronous whirl component that has limit 
cycled. They also indicate that the steady-state displace- 
ment of the rotor due to the hydrodynamic sideload is much 
larger at the turbine end of the rotor than at the turbine 
end. If only the bearing carriei' structure stiffness is 
considered, this result would not be anticipated from the 
location of the side load, i.e., comparable magni tides 
would be expected for both ends. However, the result is 
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explained by the additional stiffness provided by the seal 
coefficients of Table 5. 

Figures 8 (g) and (h) illustrate the motion of the 
pump-end squeeze-film damper cylinder, while figures 8(i) 
and (j) illustrate the motion for the turbine-end squeeze- 
film damper. The results of these figures are basically 
similar to the rotor motion at these locations iii figures 
8 (e) - (f ) . The large magnitude of the steady-state dis- 
placements in figures 8(i) and (j) precludes the use of a 
clearance smaller than 6^ = .003 in. Figures 8 Cg) - (j) 
demonstrate that the pump and turbine squeeze-film dampers 
are operating at eccentricities of 0.3 and 0.57, respec- 
tively. 

Figures 8 (k) and (1) illustrate the force magnitudes 
developed at the pump- and turbine-end squeeze-film dam- 
pers, respectively. The force magnitude at the pump-end 
squeeze-film damper is larger and more smoothly defined 
than for the turbine-end squeeze-film damper. This is the 
result of (a) a smaller clearance and eccentricity ratio 
for the pump-end squeeze-film damper, and (b) the synchro- 
nous motion of the pump-end of the rotor as opposed to the 
combined synchronous-subsynchronous motion of the turbine 
end. 

Figure 8 (m) illustrates the rotor displacement magni- 
tudes as a function of axial position, which confirms 
the previous statement that the motion at the pump-end of 
the rotor is smaller than the turbine-end. It also demon- 
strates the influence of the hydrodynamic side load in 
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causing large deflections at the approximate midspan of the 
rotor . 

The maximum pump and bearing reactions accompanying 
the results of figure 8 were 240 and 440 lbs, respectively. 

The required squeeze-film damper clearances for the 
"measured” bearing support stiffness results of figure 
4(b) (hgp = = 1.4 X 10^ Ib/in? = ?2 = .01) are 

fi (pump) = ,002 in 

P (45) 

fl^(turbine^ = .0035 in 

By comparison to Eq. (44) , one observes that the turbine- 
end clearance has been increased from ,003 in 'to .0035 in, 
V7hile the pump-end clearance is unchanged. The frames of 
figure 9 indicate the results of a simulation run with 
these squeeze-film damper clearances. 

Figures 9 (a) , and (b) illustrate, respectively, the 
modal coordinates ^X2* motion of these variables 

is similar to that of figures 8(a) and (b) ; howevex, 

now has a slight sxaDsynchronous component, and the subsyn- 
chronous component of qj^ 2 is more pronounced. The modal 
coordinate q^^j^ continues to be an order of magnitude 
greater than 

As in figure 8 , the motion of the rotor and the 
squeeze-film damper cylinder at the same axial position 
are very similar, hence only the motion of the squeeze- 
film dcimper cylinders is illustrated in figure 9. The 
motion of the pump-end squeeze-fiDm damper cylinder is il- 
lustrated in figure 9 (c) and 9 (d) , and continues to be 
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figure 9c 

Purtp-end squeeze-film danper cylinder motion in t±ie X~Z plane for 
= 1.4 X 10 ^ Ib/in; = Zz = .01« 

o 



figure 9d 

Punp-end squeeze-film daT^.>er cylinder orbit in the X-i. plane for 
k = 1.4 X 10^ Ib/in; = Cz " *01. 
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figure 9e 

Turi^ine-end squeeze-film dairper cylinder motion in the X-Z plane for 
k = 1.4 X 10^ Ib/in; Ci == Cz - 'Ol* 



figure 9f 

Turbine-end squeeze-film danper cylinder orbit for k = 1.4 x 10^ Ib/in 


= C2 = . 


Vdjia Uf 1*4 

Tlht iSECC/li5J 


' tXiO-» J 


figure 9g 

I 

Punp-end squeeze-film danper reaction imagirltude for k = 1.4 x 10' 

Cl = Cz = *01. 
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figure 9h 

Turbine-end squeeze-film diairper reaction magnitude for k = 1.4 x 


Ib/in; 


10^ Ib/in 
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figure 9i 

Eotor displacQjnent nagnitudes versus rotor axial position for 
= 1,4 X 10^ bl/in; = ?2 = *01. 



figure 10a 

5 5 

versus time for k^ = 1.4 x 10 Ib/in? k^j, = 2.8 x 10 Ib/in; 
Cl = Cz = *01. 
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basically synchronous, although a subsynchronous component 
is evident. The motion of the turbine-end squeeze-film 
damper is illustrated in figures 9 (e) and (f ) . The sub- 
synchronous component of this motion is now comparable 
to the synchronous component. The orbit amplitude of fi- 
gure 9 (f ) is somewhat larger than the corresponding orbit 
of figure 8 ( j ) . The eccentricity ratios for the puinp- 
and turbine-end squeeze-film dampers are 0.3 and 0.71, 
respectively . 

The force amplitudes for the pump and turbine squeeze' 
film dampers are given as a function of time in figures 
9Cg) and (h) , respectively, The amplitudes of both these 
force.-’ are seen to have incx'eased in comparison to figures 
8 Ck) and (1). This result is explained by the fact that 
the rotor of figure 4 Cb) is less stable than the rotor of 
figure 4 (c) , and larger damping forces are required to 
contain it. The damper force at the turbine squeeze-film 
damper is in fact approximately twice as large in figure 
9(h) as in figure 8(1), despite an increase in squeeze- 
film damper clearance from .003 in to .0035 in. However, 
this result is explained by the associated increase in 
eccentricity ratio from 0.57 to 0.71. 

Figure 9(i) illustrates the rotor displacement ampli- 
tudes versus axial rotor position. ^\Tien compared to fig- 
ure 9 (m) , this figure explains VJhy the clearance for the 
pump squeeze-film damper could remain at .002 in, despite 
a reduction in bearing support stiffness from 2.17 x 10^ 
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Ib/in to 1,4 X 10^ Ib/in. Specifically, the amplitudes of 
the forward bearing continuas to be small because the 
pump~end of the rotor has "crossed-over" the bearing cen- 
terline, and is nov/ 180° out of phase with the turbine 
end of the rotor. The increase in rotor deflections cau- 
sed by reducing the bearing support stiffness is primarily 
felt in the turbine end of the rotor. 

One V 70 uld not necessarily be convinced from the re- 
sults of figure 9 that tl>e unstable motion of the rotor 
has in fact been contained. Plowever, the final values 
from the results of figure 9 v;ere used as initial values 
for a second simulation run, and the results v;ere basical- 
ly unchaiaged , 

An inspection of the results of figures 8 and 9 de- 
monstrates that the turbine end of the rotor is more dif- 
ficult to stabilize than the pvunp end, in part, because 
the amplitudes are larger and the turbine-end squeeze- 
film damper must accordingly have a larger clearance 
than the pump end. With this in mind, the bearing sup- 
port stiffness pair 

kgp(pump) = 1.4 X 10^ Ib/in (4 6) 

kg^ (turbine) = 2.8 x 10^ Ib/in 

was investigated, since the increased turbine-be-> cing- 
support stiffness would presumably reduce the turbine-end 
amplitudes of the rotor, and allov/ a reduction in 6^, This 
change in stiffness had a negligible direct influence on 


82 


rotor stability. An unanticipated consequence of the 
change in stiffness of Eq. (46), v/as that the pump-end 
amplitudes were increased, and the pump-end squeeze-film 
damper clearance had to be increased. The resultant clear- 
ances were 

6p(pump) = .0025 in 
6^ (turbine) = .003 in 


Figure 10 illustrates the result of a simulation run with 
these clearances, the stiffnesses of Eq. (46), and 
Si “ “ -Ol* ^he initial conditions for this run were 

obtained from a previous run for which sufficient external 
damping vi/as provided to stabilize the rotor. 

Figures 10(a) and (b) illustrate the modal coordi- 
nates f qjf2' comparison to the previous results for 

these variables in figures 8 and 9, both are seen to show 
a pronounced increase in their subsynchronous whirl compo- 
nent.' This is particularly true for the variable qj^2* 

Note that q^j continues to be an approximate order of mag- 
nitudti larger than • 

The motion for the pump-end and turbine-end squeeze- 
film damper cylinders are illustrated in figures 10(c), 

(d) and 10(e), (f ) ; respectively. The motion is clearly 
unstable, with the subsynchronous component increasing ex- 
ponentially. The damper forces for these squeeze-film 
dampers are illustrated in figures (g) and (h) , and are 


also .^een to be divergent. 

One could argue that the results of figure 10 are not 
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figure 10b 

5 5 

3 ^ versus time for Jc^ = 1*4 x 10 Ib/in; = 2.8 x 10 Ib/in; 
Ca = Ca = ‘Ol* 
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figure lOc 

Pump-end squeeze-film danr:>er cylinder motion in the X-Z plane for 
k = 1.4 X 10 ^ liyin; = 2.8 x 10^ Ib/in; = C 2 = . 01 . 
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figure lOd 

Punp-end squeeze-film danper orbit for k.^ = 1-4 x 10 Ib/in.? 
= 2.8 X 10^ ]b/iu? = Ca = -Ol' 
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figure lOe 

Turbine-end squeeze-film damper cylinder motion in the X-2. plane for 
k “ 1.4 X 10^ Ib/in; k = 2.8 x 10 Ib/in; = Ca = -Ol* 
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figure lOf 

Turbine-end squeeze-film danger cylinder orbit for k = 1.4 x 10 


k , = 2.8 X 10^ Ib/ln; Ci = ?2 = •01* 
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figure lOq 

I 

Pun^end squeezer-film datnoer reaction magniti:de for k = 1.4 x 10' 

- ^ 

k , = 2.8 X 10^ Ib/in; Ci = C 2 = • 


Ib/in; 


Ib/in 
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figure lOh 

Turbiner-eud squeeze-film dairper reaction magnitude for k 

5 ^ 

= 2.8 X 10 Ib/in? Zi = = *01. 


= 1.4 X 10 lb/3Ji; 
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figure 11a 

A cxantinuation of figure lOe. Ikirbine-end squeeze-film dan^per cylinder 
motion in the X-Z plane for k„ = 1.4 x 10^ Ib/in; k . = 2.8 x 10"^ Ib/in; 
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conclusive concerning rotor stability/ since the rotor am- 
plitudes could simply be approaching a large , but bounded- 
amplitude limit cycle. To obtain more conclusive results/ 
final values from the run which produced the results of 
figure 10 were used as initial values for a second run. 

In this second run, the amplitudes continued to increase 
until the ptimp-end squeeze— film damper clearance v/as ex- 
ceeded. Figures 11(a) and (b) illustrate the motion of 
the turbine-end squeeze-film damper cylinder from this si- 
mulation run. The unstable subsynchronous whirl in figure 
11(a) is seen to be at 203 Hz. (12,200 rpm, or 1275 rad/ 
sec). Hence, the unstable motion is clearly associated 
with the first mode, and the second mode is apparently 
stable , 

The proposed squeeze-film damper redesign of figure 7 
entails a reduction in flow area in annular flow passages 
of the pump, which has the undesirable side effect of con- 
stricting coolant flow. To circumvent this problem, a 
"partial" squeeze-film damper design v/as considered which 
incorporated small radial clearances over the half-circum- 
ference of the damper which had reduced steady-state clear- 
ances due to the hydrodynamic side load and large (.008 in) 
radial clearances over the portion which had increased ra- 
dial clearances. This approach was follov’ed for the tur- 
bine-end damper for each of the "fixed-clearance" situations 
of figures 8, 9, and 10. Specifically, simulation runs v/ere 
performed with the large-clearance half circumference 


of the damper increased to .008 in, as compared to .003 in 
or .0035 in for the reduced-clearance half circumference of 
the damper. This "partial" squeeze-film damper design ap- 
proach had no discernible effect on rotor stability. It 
does, however, have the practical advantages of increasing 
the supply of fluid, for the damper , and allov/ing free cool- 
ant flow through the damper. 

A review of the results of figures 8, 9, and 10 indi- 
cates that the squeeze-film damper designs of figure 7 
can be of some value in controlling a rotor instability; 
however, the results are neither particularly encouraging 
nor conclusive. The efficacy of the dcimper is very de- 
pendent on the proper choice of clearance. Without some 
measured data to determine the amplitude of bearing car- 
rier motion, it is going to be difficult to do more than 
guess at appropriate clearance values. Figures 10 and 11 
illustrate the result of selecting damper clearances 
which are too large. Conversely y the selection of clear- 
ances which are too small would cause the dampers to 
"seize" or "lock-up", eliminating the soft-bearing design 
with its desirable critical-speed locations. 

The simulation runs v/hich yields the results of fi- 
gures 8, 9, 10, and 11 were made on an IBM 360-165 compu-^ 
ter. Each run involved 1010 integration steps with a con*^- 
stant integration stepsize of 5 x 10“^ seconds, and yiel- 
ded 17 calcomp plots in addition to periodic printed out- 
put, The computer time requirements for each run v;as ap- 
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proximately 1 minute for execution and 16 seconds for com" 
pilation. 

4.5 TRriNSIENT SIMOLATIOcJ OF HPFTP MOTION WITH SQUEEZE-FILM 
DJriMPERS AND ACYMMETRIC SUPPORT STIFFNESS. 

The results of Chapter III show that stiffness asymme- 
try in the bearing supports can improve rotor stability 
characteristics, and the preceeding section shows that 
squeeze-film damper forces can also improve rotor stability. 
The efficacy of a combination of these two approaches is 
considered in this section. The stiffness asymmetry set of 
Eq. (36) was considered for the follo\i7ing radial damper 
clearance choices : 

(a) Pump-end damper - 

reduced clearance half circumference = .002 in 
enlarged clearance half cir'cumference = .008 in 

(b) Turbine-end damper . 

reduced clearance half circumference = .003 in 
enlarged clearance half circumference = .008 in 

The simulation model was initially run with the internal 
damping factors ~ sufficient external dam- 

ping to ensure stability, and the final conditions of this 
run were then used as initial conditions for "stability" 
runs. The first such stability run was inconclusive, in 
that the rotor amplitudes increased, bvit it was not clear 
whether they would diverge or approach a stable limit. A 
subsequent run v/as made using the final values of the se- 
cond run as initial conditions. The results of this run 
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are contained in figure 12 , and clearly indicate rotor in- 
stability. Figures 12(a) and (b) illustrate the modal co- 
ordinates and ^ subsynclironous component is evi- 

dent in both these figures. Motion in the X-Z plane for 
pump and turbine bearings are illustrated in figures 12 (c) 
and (d) . The turbine-bearing motion clearly shov;s the un- 
stable nature of the motion with a whirl frequency at ap- 
proximately 12,000 rpm. Hence, it is (again) the first 
mode that is unstable. 

Although not illustrated, the simulation run which 
yielded the results of figure 12 were repeated with ? i = 

5 2 “ 0,0. For these damping coefficients, the motion is 
stable, although it cc itains a substantial subsynchronous 
component. 

The results of this section indicate that the effec- 
tiveness of a redesign which combines the stiffness asyiiuue- 
try of Eq. (36) and the squeeze-film damper design of fig- 
ure 7 depends on the degree of rotor instability. Specifi- 
cally, the addition of squeeze-film dampers can control the 
predicted instability of curve (a) in figure 4(a), but is not 
sufficient enough to control the predicted instability of 
curve (b) . 

4,6 HPOTP SIMULATION DEVELOPMENTS 

The HPOTP transient simulation model of Ref. [1] was 
extended in this study to account for the seal coefficients 
of Tables 13-16. Also, the two-mode model of Ref. [1] was 
replaced here v^ith a model for which four modes are re- 
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figuL*e 12 c 

Punp-end squeeze-film damper cylinder motion in the X-Z plane for the 
support stiffness definition of Eq, (36) vdth Ci - ?2 *= *01. 
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figure 12 d 

Turbimrr-end squeeze-film damper cylinder motion in the X-Z plane for the 
support stiffness definition of Eq. (36), vd-th = C 2 = *01. 
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tained for integration/ and an additional 9 modes ‘ sta- 
tic contribution are accounted for via a residual flexi- 
bility matrix. This improved model has been developed 
and verified; however/ no specific use has been made of 
it in the present study. It is available for use in re- 
solving rotordynamic problems v/hich the SSME test program 
may uncover. 

If the HPOTP should have a subsyhchronous whirl pro- 
blem as suggested by figures 5 the only design modifica- 
tions which this investigator could presently propose 
for study would be modifications of the clearances be- 
tv/een either (a) the outer races of the ball bearings and 
the bearing cartridge/ or (b) the bearing cartridge and 
the pump housing. It is possible that the clearances at 
these locations cguld be used to develop stabilising 
squeeze-film damper forces v/ith liquid oxygen as the wor- 
king fluid. Vance [18] discusses a squeeze-film damper 
design formed at the outer race of a ball bearing; hovj- 
ever/ his damper design is similar to that of figure 7 in 
that both an external oil supply and axial-flow seal con-, 
straints are required. No damper centering spring is 
attached to the ball-bearing outex'-race in this design, 
hov/ever. Rocketdyne personnel have suggested that the 
present clearances will reduce the rotor amplitudes during 
critical speed transition below those predicted by this 
investigator [2-3] . Test results should be available 
shortly to settle this question and to determine whether 


a closer study of the potential value of these clearances 
is merited. 

The model cited in the first paragraph of this sec- 
tion is not appropriate for the investigation of possible 
squeeze-film damp forces at the bearing outer-race clear- 
ances or bearing-cartridge clearances . Two additional 
HPOTP transient simulation models have accordingly been 
developed, based on free- free rotor modes, v/hich can be 
used Cif needed) to study possible sgueeze-fi3.m forces as 
sociated with these clearances. 
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CHAPTER V 

CONCLUSIONS AND RECOMMENDAT EONS 

Before summarizing the conclusions of this study, 
certain (possibly evasive) disclaimers are in order. The 
results of this study are at present generally without 
direct experimental verification. The study is almos!: 
entirely analytic in nature, and has entailed many assump- 
tions and many points of judgement. An effort has been 
consistently made to specify the assumptions involved, 
and to provide sufficient data to allow interested par- 
ties to check the results , and to possibly x*each conclu- 
sions different from those stated. The raw data used in 
this study have been supplied by Rocketdyne and NASA per- 
sonnel based on best-available pre-test estimates. There 
has been a concious effort made throughout this study to 
avoid reaching dire conclusions based on a series of se- 
quential ''worst-case" assumptions. The study has been 
conducted to (hopefully) anticipate rotordynamic problems 
prior to testing, and to examine some alternative means 
for rectifying these problems should they occur. It has 
had the coroll£iry objective of providing an improved un- 
derstanding of the turbopump x*otordynamic characteristics , 
v/hich should increase the probability of a correct diag- 
nosis of rotordynamic problems from test data. 

The results of this study indicate a serious rotor- 
dynamic instability problem for the HPFTP . The predicted 
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instability is caused by hydrodynamic seal forces and in- 
ternal rotor damping. 

A redesign of the bearing carrier structure to yield 
support stiffness asymmetry was considered, and substanti- 
ally improves the stability characteristics of the HPPTP 
rotor by elevating the predicted onset speed of instability. 
However, an onset speed of instability continues to be pre- 
d:t cted within the operating range for the degree of asymme- 
try which was used and assumed to be feasible in this study. 

A redesign of the bearing carrier structures to yield 
squeeze-film damper surfaces (see figure 7) was investiga- 
ted with mixed results. These dampers do generally improve 
rotor stability, but their sufficiency in this regard de- 
pends on the degree of rotor instability, and they can be 
overloaded . 

A combination redesign incorporating support- stiffness 
asymmetry and squeeze-film dampers was also examined, again 
with mixed results. Specifically, a combination redesign 
is sufficient to control a degree of rotor instability, but 
can be overloaded. 

In the event that a rotordynamic instability is encoun- 
tered in the HPPTP test program, v;hich is similar to that 
predicted in this study, the following recommendations are 
made ; 

(a) The bearing carriers should be redesigned to obtain 

stiffness asymmetry, and testing should then be conduc- 
ted to determine the efficacy of this approach. 
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(b) Only if recommendation (a) is inadequate to control the 
instability should the proposed squeeze-film damper re- 
design of figure 7 be considered. A decision to imple- 
ment the proposed redesign of figure 7 would necessar- 
ily need to be based on high quality test data, prefer- 
■ ably proximity gau-'e data of actual rotor motion. 

The results of this study concerning the HPOTP are not 
conclusive. A rotor instability associated with the first 
mode is possible, but is significantly less likely than for 


the HPFTP. 
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APPENDIX A; MODAL ROTORDYNAJIIC FORMULATIONS 

Childs in Ref. [1] developed the following two flexible ro- 
tor formulations based on vhe traditional Jef fcott-Green 
flexible rotor representation; 

{a) a modal formulation based on the "conventional" Jeff- 
cott flexible rotor model in V7hich rotor deflections 
and rotations are defined in a stationary coordinate 
system, and 

(b) a modal formulation for v/hich the rotor deflections 
and rotations are initially defined in a rotor-fixed 
cooj.'dinate system. 

In the present study, the "rotor- fixed" model was used for 
the transient simulation v?ork, while the "conventional" mo- 
del was used in stability analysis. A summary of the es- 
sential elements of these models follov/s. 

V 

A "CONVENTIONAL" MODAL JEFPCOTT MODEL 

The rotor is modeled as a collection of n rigid bodies con- 
nected to each other by an elastic structure. The rotations 
and deflections of component rigid bodies are defined with 
respect to an inertial X, Y, Z coordinate system, v^here Z 
defines the rotor's nominal axis of rotation. The original 
structural-dynamic model for the rotor can be put into the 


form 


102 


tm] o' 
0 [J] 

,(0y)J 

[•'• '"SezJ ^ 

(&yU 

Im] o' 
0 [J]_ 


j.+ < 

f<v', 

i«X>J 






(*y) 

‘”bc> 


(A.l) 


L + 


v/here (R^) , (Sy) and (R^) , (g^^) define the deflections and 

rotations in the X-Z and Y-Z planes, respectively. Further, 
[m] and [J] are the diagonal mass and diametral-moment-of- 
inertia matrices. The stiffness matrices [1^2^ ' de- 

fine the nominal axisyiranetric' linear stiffness properties 
of the rotor-bearing structure. The vectors (fj^) , (m^) and 
(f^) f (m^j) are the external forces and moments in the X-Z 
and Y-Z planes , respectively. 

Conventional eigenanalysis of the left hand side of equa- 
tion (Ail) yields 2n eigenvalues ^ and 2n eigenvectors 
(A^) . The eigenvectors define the follov/ing coordinate 
trans formations 

(%) = [A^lCqY)/ (3x^ =-[Ap]tgY> 


(A. 2) 


From equation (A. 2), tlae complete eigenvector matrices for 
the X-Z and Y-Z planes are defined by 


[Ajjg] 


Ao 




e 

J 


The eigenvector matrices are normalized to satisfy 


(A. 3) 


[Ag] + tAp]^[J] [Ag] = [D] 

IAxz]^[Kx2^ 

vrhere [U] is the unity matrix^ and [A] is the diagonal matrix 
of eigenvalues . 

The coordinate transformation defined by equations (A. 2) 
and (A. 3) yields the following modal vibration equations in 
the X-Z and Y-Z planes 

(gjj) = " ^[CM](gY) - [AHg^) t *>■ i tPy) - 

[C]{(4^)+ ^(qy)) 

(A. 4) 

(q^) = (Fy) + I ICM] (q„) - [ A] (qy) + J (p^) _ 

[C] { (q^)“ ij) ^^X^ ^ 

v/here <|) defines the rotation of the rotor about the Z axis , 
and 

[CM] = [Apl^’tJ] [Ap] 

(P^) = (Px) - s4> (Py) , (Py) = (P^) + C((. (Py) 

(Pj.) = [Agl'^Cma^) + (A,5) 

(Py) = [A^]^(raay) + [Ap]^(Jy2) 

In the above, [j”] is a diagonal matrix of polar moments of 
inertia and the components of the vectors (ma^) , (mSy) , 

' CJy 2 ^ ^respectively , where 

a^j,, a^y define the unbalance of component rigid body i, and 
s^ire the products of inertia of rigid body i. 
Returning to equations (A. 4) , the vectors (Pv ) r (Pv^ 
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are generalized modal forces defined by 

= tAg]®(fjj) + [Ap]^(My) 

(Py) = [Ag]^(fY) " 


(A. 6) 


The matrix [C] in Eg. (A. 4) defines the internal damping of 
the rotor, and has the diagonal elements where is 

the damping factor. Component equations' from (A. 4) are 


P;gjL “• <f3:CMj^jgyj - 



(qjji + «Jqy£) 

-1- t iPy^; 

i = 1 , 

2, , . 

. 2n 




<A 

Pyi t I^CMijg^. - 

Xi^qyi 

- 2?i\ 

(4yi - iqxi> 


i == 1, 

2 ^ « * 

. 2n 


The governing equation for (|» is 
jJ = (Py)^Cqx^ 


(A. 8) 


where J is the polar moment of inertia for the rotor, and 
M is the resultant spin-axis torque applied to the rotor. 

A procedure for accommodating nonlinear or asymitietric 
bearings in this analysis can be explained by the following 
example. Suppose the bearing reaction at the ith station 
is defined by the components 

q/ - 

Ry^, Ryi, r/, <f) 

A Taylor series expansion of these functions about R^jj = 

= 0 , I = 4» yields a reaction definition of the form 


(A. 9) 


>x"’ 




- = ^ 


+ - 





gx^ (3^^ , Ry^ , , Ry^ r 5 I ) 

[gy^ (Rjji , (Ry^ f Rjji , Ryi , 6 ) 


{A. 10) 
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The matrix [K^] can be stated 
[K^] = [K^^l + tK^^g] 


(k 1 1 + kg2 ) 0 

0 (3^22 + k 1 1) /2 

(kii- k22)/2 0 

0 (k2 2 ” J^ii)/2 


(A. 11) 


The matrix [K^_] is absorbed in the system stiffness matri- 
a 

ces [Kjjg] f of equation C8) , and the reaction terms 

of equation (20) become 


q/ 


- [Kg^] 




’ i 


i i 

. > 

a ^ 


(A- 12) 


The negative of these terms are treated as external forces 
acting on body i^ and are accounted for in the generalized 
force vectors of equation (A, 8). Note that this develop- 
ment does not linearize the bearing characteristics. Spe- 


f account for the asymmetric, 


cifically, the terms 
nonlinear , and speed dependent properties of the beajrings . 
A similar procedure is followed to account for forces ari- 
sing from seals and fluid-film dampers. 


A "Rotor-Fixed*' Modal Jeffcott Model 

The equations of motion to be presented here differ from 
the preceeding in that they are stated in a rotor-fixed x, 
y, z coordinate system. The modal equations are 

(q^) = (F^) + <J{2[U] - [CM]}(qy) - {[A] + <j: 2 [CK] } (q^) 


X 

[C] (q^) + (P^) -!■ 


(A. 13) 


(q^) = (Fy) - <^{2[U] - [CH]}{q^) - { [ Al + [CK] } (q^) 
** [C] (q^) - J(P^) + 


? 1 
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vj-here 

(Fy) = lAgl'^Cfy) - lAp]"^ (A,14) 

[CK] = [CM] ~ [Ul 


The component fom of equation {A. 13) is 


2t,.l .q . + 4P . + A^p y i, = If 2, . 
X i^xx ^ yx XX 


. 2n 


(A. 15) 


-Jyi = ^yi - - ™ij)9x3 - 


2i?x.a . ~ AP . + <j2p .. i = 1, 2, 

■i x^yx ^ XX ^ yx '^ , 


.2n 


where 5 . . is the Kronecher delta* 

The transformation from the rotor-fixed modal coordinates 
of Eq. (A. 15) to the displacements and velocities of compo- 
nent rigid bodies is 


(Rx^ = [AgHc<i.(q^) - S(j.{qy)} 

(Ry) = [Ag]{s<^(q^) + c<|.(gy)} 

(A. 16) 

(^) = [A^] {c<f. (q^ - .^qy) - S(f> (qy + jq^) } 

C^) = [Ag] {S(j) (q^ - <fqy) -f ci> (qy + (jiq^) } 
and the differential equation for 4i is 

jJ = - 2<j;{ (P^)^(4^) + (P^)^(qy)} + {Py)^{q^) 

- " <^^UPy)^(q^) - CP^)^(qy)> (A, 17) 

The procedure followed in past analyses has been to 
retain for integration those modes whose natural frequen- 
cies are less than the rotor’s top running si^eed plus one 
or two modes above the top running speed. However in the 
present study, the residual flexibility method of Schwend- 
ler and Macneal was adapted to account for the ’’static" 
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contribution of higher order inodes. The technique employed 
is simple r and can be expla.uied by considering the following 
structural-dynaiviic modal 
[M] (X) + [K] (X) = (f) 

where [M] and [K] are n x n mass and stiffness matrices, 
and (X) and (f) are n x 1 vectors of generalized forces and 
displacements. Eigenalysis of this system yields a matrix 
of eigenvectors [A] , and a diagonal matrix of eigenvalues 
[ A] v/hich satisfy 

[A]^[M][A3 = Etj], [A]^[K][AJ = [A] 

Further, the coordinate transformation 

(X) = EA] (q) (Avl8) 

yields the modal differential eqxxations 

(q) + lA](q) = tA]’'(f) (A.19) 

The first n modes are to be retained for integration, and 
the static contribution of the x-emaining (n-m) modes are 
defined by the algebraic equations 

[ A] (g) = [A]^{f) (A. 20) 

where the elements of I A] are the (m + 1) through n eigen- 

m 

values, and the matrix [A] contains the (m + 1) through 
n eigenvectors. Denoting the m modes to be retained for 
integration by q yields the following complete model 

(t) t E:a] (q) = [A]^(i) (A. 21) 

(x) = EA]^{q) + EA) (q) 

= EA]^(q) + EA] EA3''^EA]^{f) 

= EA]^(q) + E^Hf) 

The matrix E?’3 is called tlie residual ^flexibility matrix and 
accounts for the static contribution of higher order modes. 
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APPENDIX B SSME TURBOPUMP ROTOR/BE?VRING DYNAMIC MODELS 

The structural dynamic models for the SSME turbo- 
pumps v/ere developed by Rocketdyne personnel, and are 
discussed in detail in References {!] , [2] , with geome- 
tric, inertial, and eigen data provided. Hence, only a 
brief review of these models is given here, and only 
eifferent or additional data is provided. 

HPPTP MODEL 

The HPFTP jotor is modeled by twenty seven eLasti- 
cally connected rigid bodies while the bearing outer 
races and carrier structures are modeled by two point 
masses. In the model, the bearing support mass is at- 
tached to '-'ground" via a linear spring with coefficient 
kg, and is attached to the rotor through a ball bearing 
of stiffness kj^ as illustrated in figure B.l, The de- 
sign value for kg is 

kg ^ 2.17 X 1Q5 Ib/in (B.l) 

The same bearing is used for both the pump- and turbine- 
end bearings, and the stiffness of this bearing is defined 
as a function of running speed | by 
kj^ - Uq + a^l + 
ao = 1.475 X 106 Ib/in 

aj = 162,8 Ib/in (B.2) 

a 2 = -.1559 lb sec^/in 
as = 1,949 X 1Q6 Ib/in 
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The basic eigendata used in this study was calculated for 
the support stiffness of Bq. (B,l) , and a nominal bearing 
stiffness = .891 x 10 ® Ib/in, which is the approximate 
value of kj^ in Eg, (B,2) at FPL, The first four eigenvec- 
tors are given in Table B.l, and correspond to the fol- 
lowing natural frequencies . 

Al = 921.5 rad/sec = 8,800 rpm 

^2 = 1589. rad/sec = 15,200 rpm 

A 3 = 3558. rad/sec = 34,000 rpm 

A 4 = 7486. rad/sec = 71,500 rpm 

In addition, the next eight higher eigenvalues and 
their associated eigenvectors were used to calculate a 
residual flexibility matrix follovjing the procedure outlined 
in Apx>endix A. The residual flexibility matrix was used 
in tlie sijiiulation model, but v?as not used in stability stu- 
dies. The residual correction technique v?as required to 
obtain the correct static load-def 3 ection characteristics 
for the bearing support structure. Specifically, the first 
four modes of Table B.l are not adequate to define the cor- 
rect static load-deflection characteristics for the coor- 
dinates ^Sl' ^b2'^ ^S2 figure B.l. The first four 

modes give reasonably good results for the displacements 
along the rotor and at the bearings in 

error for both the bearing-support mass coordinates 
Xg 2 and the relative motion between the shaft and the 
bearing-support masses. The explanation for this situa- 
tion is that the eigenvectors which define the bearing sup- 
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port mass degree of freedom relative to the rotor are as- 
sociated with natural frequencies that are much higher than 
the first four modes of Eg. (B,3). Correct static load- 
deflection properties v;ere obtainable from the first eight 
modes. However (as noted above) the first tv/elve modes 
were employed in the transient model with the first four 
modes retained for integration and the next eight modes 
used to calculate the residual f lexibility ^matrix . 

HPOTP MODEL 

The HPOTP rotor is modeled by thirteen elastically- 
connected rigid bodies. The HPOTP bearings are supjported 
directly by the turbopuxnp housing in contrast to the soft- 
bearing-mount design used for the HPPTP. The bearing 
stiffness definitions used for the pump and turbine 
bearings are 

kp(pump) = a^Q + t + a^^l^ 

SpQ = 1.321 X 10^ Ib/in 
a^T - 81.28 lb sec/in 

(B.4) 

3p2 “ -.1857 lb sec^/in 

a ^ = 3,633 X 10*"^lb sec^/in 
P3 

(turbine) = a^^ + ®t3^^ 

a^Q = 1.79S X IQG Ib/in 
a^j^ = 105.2 lb sec/ in 
a ^2 ~ -.3124 lb sec^/in 
a ^2 " 5.912 X 10“^ lb sec^/in 


(B.5) 
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Eigenanalysis v/as performed for the nominal bearing stiff- 
nesses k = jc . = .891 X 10® Ib/in- The first four eigen- 
p t 

values for these nominal bearing stiffnesses are 
Xj = 1201. rad/sec = 11,463 rpm 

A2 = 4082. rad/sec = 38,985 xrpm 

A 3 = 5946. rad/sec 56,780 rpm 

Aif = 6440. rad/sec 61,500 rpm 

Their associated eigenvectors are given in Table B.2. 

Only the first four modes were used in the stability an- 
alysis of the HPOTP; however, the next nine modes v?ere 
used to calculate a residual flexibility matrix for the 
transient simulation model. 


Zi 


[Ag] X 10 


S O 

*=51 W 

sg 

H ^ 

^63 


tA^l 


1 . .9960 

0.9118 

“ 3.191 

- 5.264 

- 6.151 

- 1.077 

- 1.579 

“ 7.339 

- 15.67 

2 . 2 . 29 S 

1,0522 

- 2.986 

- 4-303 

- 4.031 

- 1.067 

- 1.573 

“ 7.046 

- 13.89 

3 . 2.946 

1.122 

- 2.884 

- 3.825 

- 3.083 

- l.OSO 

- 1.576 

“ 6.932 

- 13.11 

4 . 3.646 

1-198 

“ 2.774 

- 3.318 

- 2.053 

- 1.049 

- 1.578 

“ 6.721 

- 11,85 

5 . 4.516 

1.302 

- 2.662 

- 2,737 

- 0.9417 

- 0.9504 

- 1 . S 57 

- 6.322 

- 9.230 

6 . 6.416 

1.492 

“ 2.353 

“ 1.560 

0.6904 

“ 0.9153 

“ 1.781 

“ 6.057 

- 7.187 

7 . 7.966 

1.635 

- 2.076 

- 0.6028 

1.645 

“ 0.8121 

“ 1.923 

“ 5.423 

- 3.919 

8 . 9.096 

1.728 

- 1.856 

0.0169 

1.944 

- 0.7007 

- 2.064 

- 4.172 

- 1.005 

9 . 11.40 

1,871 

- 1.355 

0.9958 

1.550 

- 0.5426 

- 2.247 

- 3.664 

2.413 

10 . 12,95 

1.949 

“ 0.9946 

1.513 

0.8412 

- 0.4231 

“ 2.333 

- 2.627 

3.424 

11 . 14.16 

1.985 

“ 0.6963 

1 . 719 - 

0.1253 

- 0.1871 

“ 2.487 

- 0.7536 

4.388 

12 . 16.38 

1.996 

“ 0.1202 

1.631 

“ 1.135 

“ 0.0177 

“ 2.587 

0.9279 

4.064 

13 . 18.05 

1.978 

0.3211 

1.341 

“ 1.766 

0.2111 

- 2.596 

2.427 

0.9768 

14 . 19.57 

1.934 

0.7201 

0.8786 

“ 1.810 

0.5635 

- 2.569 

3.360 

- 1.559 

15 . 21.00 

1.871 

1.094 

0.2740 

“ 1.415 

0.5116 

- 2.472 

4.036 

- 3.871 

16 . 22.25 

1.807 

• 1.414 

- 0.2870 

- 0.8413 

0.4993 

- 2.532 

4.346 

- 4.027 

17 . 23.50 

1.738 

1.727 

- 0.363 

- 0.1280 

0.5178 

- 2.519 

4.410 

- 5.229 

18 . 24.88 

1.658 

2.070 

- 1.508 

0.7132 

0.5317 

“ 2.507 

4.458 

- 5.512 

19 . 26.05 

1.584 

2.345 

- 2.022 

1.388 

0-5817 

“ 2.436 

4.446 

“ 5,704 

20 . 26.62 

1.539 

2.467 

- 2.268 

1.737 

0^6048 

- 2.405 

4.456 

“ 5,856 

21 . 27.63 

1,476 

2.708 

- 2.726 

2.370 

0.61141 

“ 2.399 

4.474 

- 5.969 

22 . 4.246 

1.260 

“ 2.690 

- 3.020 

- 1.666 

- 1.067 

“ 1.512 

“ 6.497 

- 1,174 

23 . 7.496 

1.591 

- 2-163 

“ 0.9348 

1.350 

- 0.9234 

- 1.806 

“ 6.211 

“ 5.932 

24 . 10.40 

1,822 

“ 1.588 

0.6782 

2.053 

- 0.6562 

- 2.150 

- 4.560 

1.007 

25 . 12.88 

1.949 

- 1.017 

. 1.528 

1.046 

- 0.3601 

- 2.439 

- 2.184 

6.396 

26 . 15.35 

1.994 

- 0.4035 

1.784 

- 0.7409 

- 0.0819 

- 2.601 

0.2458 

6.487 

27 . 17.75 

2.005 

0.2402 

1.386 

- 1.830 

0.1976 

- 2.642 

2.429 

2,023 

28 . 3.646 

0.9611 

- 2.236 

- 2.751 

- 1,952 





29 . 26.62 

1.234 

1.936 

- 1.867 

1.581 






TABLE B.l HPFTP Rotating Assembly Eigenvectors. 
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Zi [Ag] 


1. 1.875 

.0530 

■1.030 

7.372 

3.910 

2. 2.920 

-.0304 

■.2070 

5.669 

2.475 

3. 3.930 

-.1110 

.5361 

3.985’ 

1.075 

4. 4.800 

-.1847 

1.275 

2.879 

.0386 

5. 6.430 

-.3187 

2.469 

1.875 

-1.158 

6. 9.600 

-.3337 

3.352 

-9114 

“1-059 

7. 13.20 

-.2052 

2.982 

-.4729 

.3777 

8. 14-75 

.1795 

2.301 

-.7391 

1.178 

9. 15.87 

.5329 

1.974 

-.8999 

1.720 

10. 16.93 

.9266 

1.741 

“1.027 

2.177 

11. 18.37 

1.527 

1-448 

“1.093 

2.547 

12. 20.10 

2.349 

1.060 

-.9141 

2.313 

13. 22.25 

3.404 

.4901 

-.2110 

.C529 

TABLE 

B.2 First 

Four 

Eigenvectors 

HPOTP^ 


[A,] 


.0798 

-.7844 

.0798 

-.7818 

.0797 

-.7684 

.0779 

-.7368 

.0337 

-.3690 

.0074 

-.1227 

.1974 

.2498 

.2597 

.2899 

.3400 

.2622 

.3794 

.2399 

.4244 

.2337 

.4704 

.2648 

.4932 

-3131 


1.569 1,338 

1.549 1.322 

1.392 1.218 

1.133 1.032 

.4645 .1647 

.3981 “.1951 

.2140 “.4429 

.1428 “.4321 

.1073 -.3729 

.0578 -.2605 

-.0426 -.0040 

-.2367 ,5294 

-.4137 1.034 
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